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Abstract 

A  general  mixed  Hj/H^  optimal  control  design  with  multiple  H„  constraints  is 
developed  and  applied  to  two  systems,  one  SISO  and  the  other  MIMO.  The  SISO  design 
model  is  normal  acceleration  command  following  for  the  F-16.  This  design  constitutes  the 
validation  for  the  numerical  method,  for  which  boundaries  between  the  Hj  design  and  the 
H,  constraints  are  shown.  Hie  MIMO  design  consists  of  a  longitudinal  aircraft  plant 
(short  period  and  phugoid  modes)  with  stable  weights  on  the  Hj  and  transfer 
functions,  and  is  linear-time-invariant.  The  controller  order  is  reduced  to  that  of  the  plant 
augmented  with  the  Hj  weights  only.  The  technique  alio*  singular,  proper  (not 
necessarily  strictly  proper)  H.  constraints.  The  analytical  nature  of  the  solution  and  a 

numerical  approach  for  finding  subop timal  controllers  which  are  as  close  as  desired  to 
optimal  is  developed.  The  numerical  method  is  based  on  the  Davidon-Fletcher- Powell 
algorithm  and  uses  analytical  derivatives  and  central  differences  for  the  first  order 
necessary  conditions.  The  method  is  applied  to  a  MIMO  aircraft  longitudinal  control 
design  to  simultaneously  achieve  Nominal  Performance  at  the  output  and  Robust  Stability 
at  both  the  input  and  output  of  the  plant 


xiv 


Recently,  there  has  been  a  great  deal  of  interest  in  formulating  a  mixed  H/H-  control 
methodology  which  can  handle  bounded  spectrum  and  bounded  energy  inputs 
simultaneously.  Early  approaches  included  solving  the  problem  for  one  input/one  output, 
one  input/two  outputs  and  two  inputs/two  outputs  ([BH89J;  [KDGB90];  [DZB89]; 
[YBC90];  [MG88];  [MG90];  [KR91]).  The  general  formulation  of  the  mixed  em 
with  two  exogenous  inputs  and  two  controlled  outputs  was  first  approached  for  lu..  state 
feedback  in  [KR91].  Ridgely,  Valavani,  Dahleh  and  Stein  [RVDS92]  developed  a  solution 
for  the  general  mixed  ItyF L  problem  with  output  feedback  which  results  in  a  controller 

order  equal  to  or  greater  than  the  order  of  the  underlying  system  augmented  the  H2 
weighting  and  the  H.  weighting.  Also,  the  assumption  is  made  that  the  underlying 

problem  is  regular  and  has  no  feedforward  term.  Walker  and  Ridgely  [WR94a] 
reformulated  the  general  mixed  Hj/H^  problem  with  the  strictly  proper  and  regularity 
assumptions  relaxed  to  allow  singular,  proper  H„  constraints.  Furthermore,  Walker  and 
Ridgely  showed  strong  theoretical  results  for  controllers  selected  to  have  orders  equal  to 
or  greater  than  the  order  of  the  underlying  Hj  problem. 

Multiple  objective  optimal  control,  as  formulated  by  the  above,  allows  the  designer  to 
determine  the  tradeoff  between  noise  rejection  ( Hj )  and  some  unstructured  perturbation 
( FL  ),  which  embodies  desired  performance  and  margins  at  either  the  input  or  output  of 

the  plant  (or  some  combination).  However,  the  unstructured  perturbation  approach  to 
the  HL  problem  is  generally  conservative.  A  better  approach  is  to  exploit  the  structure  of 
the  perturbations  [DWS82],  but  the  single  constraint  in  the  mixed  H2/H„.  setup  is 
unable  to  do  this.  Doyle  [Doy82]  introduced  structured  singular  value  (|i)  synthesis  to 
design  controllers  which  are  less  conservative.  While  this  approach  handles  structured 


uncertainties,  the  current  p-syn  thesis  method  generally  results  in  a  large  controller  order. 
It  is  desired  to  develop  a  control  synthesis  method  which  can  reduce  the  controller  order 
below  that  of  p-synthesis.  One  approach  to  this  problem  is  to  consider  each  perturbation 
as  an  individual  constraint  and  solve  the  problem  using  mixed  Hj/H.  synthesis  with 
multiple  H.  constraints.  An  advantage  of  this  technique  is  that  it  allows  the  controller 
order  to  be  reduced  to  the  order  of  the  Hj  problem.  Further,  by  employing  mixed 
optimization,  one  can  design  a  controller  which  minimizes  the  effect  of  white  noise  inputs 
as  well  as  bounded  energy  inputs. 

This  diesis  will  first  develop  the  Hj/H.  problem  with  a  single  constraint.  Second, 
the  Ha/H.  problem  is  extended  to  allow  multiple  HL  constraints.  The  nature  of  the 
solution  will  be  compared  with  p-synthesis.  A  SISO  example  will  represent  the  validation 
of  this  new  technique  and  a  MIMO  example  will  address  control  design  requirements 
(Robust  Stability  and  Nominal  Performance). 
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n.  Background 


This  chapter  is  intended  to  lay  the  foundation  for  the  specific  compensator  designs  that 
will  follow.  Figure  2-1  shows  the  closed  loop  transfer  function  T,  and  it  is  desired  to 
minimize  an  appropriate  norm  on  T  due  to  varying  assumptions  about  the  exogenous  input 
signal  and  the  exogenous  output  signal  For  instance,  if  the  exogenous  signal  is  not 
known  exactly  but  is  known  to  lie  in  a  set  (p=2,o o),  then  a  reasonable  measure  for 
performance  is  one  which  looks  at  the  worst  possible  output  In  particular,  assume  that 
die  set  of  exogenous  inputs  is  given  by 

{we/,1  HI,  £  1};  p-  2,oo 


The  2-noim  is  the  energy,  and  the  «®-norm  is  the  maximum  magnitude  of  the  signal.  A 
good  measure  of  performance  is  given  by 


T 


Figure  2-1  Closed  Loop  system  T 
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which  is  the  norm  of  the  worst  possible  output  as  the  exogenous  signal  ranges  over  the 
allowable  set  The  controller  design  problem  is  given  by 


inf 


)=  inf  |T|,  _Mllctd 

“  K  (Ubiliztm  m 


Kiufailiang 

This  performance  objective  is  known  as  the  minimax  objective.  The  controller  is  designed 
to  guard  against  all  exogenous  signals  in  the  allowable  set  [DD93].  Table  2-1  shows  the 
lp^-induced  norm  for  different  exogenous  inputs  and  outputs. 


Table  2*1  Induced  norms 


input 

output 

II  wllj 

II  w  IL 

II  z  II2 

IITIL 

** 

II Z IL 

IITII2 

II  T(t)  II, 

**  not  induced  norm  exits 

For  more  information  on  signal  theory,  refer  to  [DD93]  and  other  references.  When 
uncertainties  are  in  the  system,  the  minimization  of  II  T  llp  [p=l  °°  (the  system  2-norm  is 
not  good  for  uncertainty  management)]  is  conservative,  especially  when  the  uncertainty 
model  is  highly  structured.  In  this  case,  ^.-analysis  is  a  better  tool  for  analyzing  the 
robustness  of  the  system.  Next,  we  examine  the  H„,  and  p-syn  thesis  design 
procedures.  Minimization  of  |T(t)Jj,  known  as  /,  optimization,  will  not  be  covered  in  this 
thesis. 


2.1  Hh  Optimize 


Hj  optimization,  which  parallels  the  popular  LQG  problem  in  the  optimal  output 
feedback  case,  is  based  cm  minimizing  the  2-norm  of  a  transfer  function  matrix  from  white 
noise  inputs  to  controlled  outputs  [DGKF89].  The  white  noise  input  is  assumed  to  be 
zero-mean,  unit  intensity,  and  possess  a  Gaussian  distribution.  Figure  2-2  shows  the  basic 
Hj  design  diagram  where : 
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Figure  2-2  H2  Design  Diagram 


z  is  the  controlled  output  (exogenous  output) 
w  is  a  white  Gaussian  noise  with  unity  intensity  (exogenous  input) 
u  is  the  controller  input  to  the  plant 
y  is  the  measured  plant  output 
P  includes  the  design  weights  and  the  plant 
K  is  the  controller 


The  H2  design  objective  is  to  find  an  admissible  (stabilizing)  K(s)  that  minimizes  the 
energy  of  die  controlled  output  (z),  which  is  equivalent  to  minimizing  the  two-norm  of 
T 

The  optimal  II  Tw  llj  is  represented  by  a0,  with  the  corresponding  K(s)sK(s)2apt.  Kfs)^  is 
unique  and  full  order  (the  order  of  the  nominal  plant  plus  the  order  of  the  H2  weights),  hi 
state  space,  the  plant  P  is  described  by: 


x  =  Ax  +  Bww  +  Buu 

z  =  Czx  +  D„w  +  Dbu 

y  =  Cyx  +  Dy*w  +  DynU 

The  following  assumptions  are  now  made: 

(i)  Dw  =  0 
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<ii)  ^=0 

(iii)  (A,  Bn)  stabilizable  and  ( Cy,  A)  detectable 

(iv)  DTaiDzu  full  rank ;  full  rank 


(v) 


(vi) 


[a  -  joji 

Bo 

[  c. 

Da 

NM 

3 

•n 

1 

< 

Bwl 

L  Cy 

DywJ 

has  full  column  rank  for  all  0) 


has  full  row  rank  for  all  co 


where  assumption  (i)  is  a  requirement  for  the  two-norm  of  the  transfer  function  to  be 
finite.  The  condition  on  simplifies  the  problem,  but  it  is  not  a  requirement.  For  a 
stabilizing  compensator  to  exist,  (iii)  must  be  satisfied.  Condition  (iv)  is  required  to  avoid 
singular  control  problems.  Finally,  conditions  (v)  and  (vi)  guarantee  the  existence  of 
stabilizing  solutions  to  the  algebraic  Riccati  equation  (AREs)  that  are  involved  in  the 
solution  of  the  Hj  problem.  For  a  complete  description  of  the  Hj  solution,  see  [DGKF89]. 


The  objective  of  H_  optimization  is  to  minimize  the  energy  of  a  controlled  output  to  a 
deterministic  input  signal  that  has  bounded,  but  unknown,  energy.  In  the  H^.  problem  the 
controlled  output  is  e,  and  the  exogenous  input  is  d;  therefore,  the  H„  problem  is 


inf 

Kadmuribte 


where 


l|TrfO<ol  s  sup  efoOrn)] 

The  optimal  controller  Kfs)^  yields  II T^  II..  =  y,  and  a  family  of  subop timal  controllers 
such  that  II  Tjj  ll„  <  y  can  be  defined,  where  7  >  y .  Figure  2-3  shows  the  block  diagram 
for  the  design.  In  this  case,  the  state  space  matrices  for  P  are: 
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HZH 

Figure  2-3  H.  Design  Block  Diagram 

The  following  assumptions  on  the  state  space  matrices  are  made: 

(i)  0^=0 

(ii)  (A,  Bu )  stabilizable  and  (Cy,  A)  detectable 

(iii)  full  rank;  D>rfDTyd  full  rank 

fA-j(Dl  B»"| 

(iv)  ^  ^  has  full  column  rank  for  all  co 

fA-jfOl  Bd*] 

(v)  I  has  full  row  rank  for  all  co 

L  Cy  DydJ 

Condition  (i)  is  not  required,  but  simplifies  the  problem.  For  stabilizing  solutions  to  exist, 
condition  (ii)  must  be  satisfied.  In  order  to  avoid  singular  problems,  condition  (iii)  is 
required.  Conditions  (iv)  ami  (v)  along  with  (ii)  guarantee  that  the  two  Hamiltonian 
matrices  in  the  corresponding  H2  problem  belong  to  dom(Ric).  Notice  that  there  is  not  a 
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restriction  on  D^,  because  it  does  not  make  the  closed  loop  infinity-norm  infinite.  For 
more  information  of  the  complete  solution,  refer  to  [DGKF89J. 

3  Structured  Singular  Value 

This  section  presents  a  short  synopsis  of  Packard  and  Doyle  [PD93].  The  system  is 
linear  time  invariant  with  complex  perturbations.  For  more  information  refer  to  [PD93]. 
Consider  M  €  C™.  In  the  definition  of  p(M)  there  is  an  underlying  structure  A  on  which 
everything  in  the  sequel  depends.  This  structure  may  be  defined  differently  for  each 
problem  depending  on  the  uncertainty  and  performance  objectives  of  the  problem.  This 
structure  depends  on  the  type  of  each  block,  the  number  of  blocks,  and  their  dimensions. 
These  blocks  can  be  repeated  scalar  blocks  and/or  full  blocks,  where  S  and  F  denote  the 
number  of  repeated  scalar  blocks  and  the  number  of  full  blocks,  respectively  (scalar  S: 
r,,...,rf;  full  block  F:  m,,...,mp).  Therefore,  A  is  defined  as 

A  —  {dlag[8j/r( ,...,Ss/rj 8,  gC, 

As+y  eC^.l <i  <S,1  <j< F} 

s  p 

and  /}  +  ^ntj—n  gives  consistency  among  all  dimensions.  The  norm  bounded 
m  y-1 

subsets  of  A  are  defined  as 

Ba  =  {AeA:6(d)<l} 

For  M  €  O”,  pA(M)  is  defined 

MM) - - : - - - 

min{o(A):  A  e  A,  det(I  -  MA)  =  0} 

unless  no  Ae  A  makes  I  -  MA  singular,  in  which  case  nA(M)sO.  Qearly,  pA(M)  depends 
on  the  block  structure  as  well  as  the  matrix  M.  In  general,  pA(M)  can't  be  calculated 
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exactly,  and  its  value  is  placed  between  lower  and  upper  bounds  for  certain  type  of  A 
block  structures  (scalar  blocks  or  complex  uncertainty  blocks;  see  [PD93]).  Two  special 
cases  of  pA(M)  are: 

.  if  A  =  {5/:8eC);  (S  =  l,  F  =  0,  rx  =n) 
then  pA(M)=p(M),  (the  spectral  radius  of  M) 

•  if  A  =  Cn“  ;(S  =  0,  F  =  1,  m,  =  n) 
then  pA(M)=C(M) 

For  all  but  the  two  special  cases  above,  jx  is  bounded  by 

p(M)  <  |XA(M)  <  a(M) 

These  bounds  by  themselves  may  provide  little  information  on  the  value  of  |x,  because  the 
gap  between  p  and  O  can  be  arbitrarily  large.  These  bounds  are  refined  with 
transformations  on  M  that  do  not  affect  jxa(M),  but  do  affect  p  and  O.  To  do  this,  define 

two  subsets  of  C"™ 

Q  =  {Q  e  A:  Q-Q  =  I.} 

D  =  [diag{D1,...,Ds,ds+;Im  ....jds+pI^JrD,-  e C  *  ,D,  >0,  £ R,  dJ+;  >0] 

Notice  that  for  any  Ac  A,  Qe  Q,  and  De  D 

Q*  eQ,  QAe  A,  AQeA, 

G(QA)  =  a(AQ)  =  a(A), 

DA  =  AD 

Theorem  3.8  from  [PD93]  says:  For  all  Qe  Q  and  De  D 

^a(MQ)  -  |J.a(QM)  =  |Aa(M)  =  pA(DMD'1 ) 

Therefore,  the  bounds  can  be  tightened  to 

max  P(QM)  <  max  p(am)  =  jha(m)  <  inf  c(dmd  1 ) 

QeQ  AeBA  DeD 
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The  structured  singular  value  is  a  framework  based  on  the  small  gain  theorem,  in  which 
the  robustness  of  a  system  can  be  quantified  [ABSB92].  n- based  methods  have  been 
useful  for  analyzing  the  performance  and  robustness  properties  of  linear  feedback  syste  m, 
where  the  closed  loop  system  and  weighting  functions  are  contained  in  the  M(s)  matrix, 
and  all  uncertainty  blocks  are  put  into  a  block  diagonal  A(s)  matrix.  M(s)  and  A(s)  are 
stable  transfer  functions;  they  are  arranged  as  shown  in  Figure  2-4.  This  figure  is  meant  to 
represent  the  loop  equations  e  =  Md,  d=  Ae.  Assuming  a  fixed  s=yco,  as  long  as  I  -  MA  is 
nonsingular,  the  only  solutions  e,  d  to  the  loop  equations  are  e  =  d  =  0.  However,  if 
I  -  MA  is  singular,  there  are  infinitely  many  solutions  to  the  equations,  and  the  norms  Hell 
and  lldll  of  solutions  can  be  arbitrarily  large;  therefore,  this  constant  matrix  feedback 
system  is  "unstable".  Likewise,  the  term  "stable”  will  describe  the  situation  where  the  only 
solutions  are  identically  zero.  In  this  context,  |iA(M)  provides  a  measure  of  the  smallest 

structured  A  that  causes  instability  of  the  constant  matrix  feedback  loop.  The  norm  of  this 
destabilizing  A  is  exactly  l/fiA(M)  [PD93].  This  interpretation  can  then  be  repeated  for  all 

frequencies. 


Figure  2-4  M-A  feedback  connection 
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The  robustness  of  a  closed  loop  system  can  be  analyzed  by  forming  the  block  diagram 
as  shown  in  Figure  2-5,  where  dA  and  eA  are  the  inputs  and  outputs  related  with 

uncertainty  block  Aj.  The  inputs  and  outputs  related  to  the  performance  specification  are 
given  by  d  and  e. 


Figure  2-5  Robust  Performance  Diagram 
The  transfer  functions  between  inputs  and  outputs  are: 


Te  ]  [Mu  M,2T  d  I 

LeJ  L^2i  M22J|_dAJ 

dA  ~  ^2e& 

This  set  of  equations  is  called  well  posed  if  for  any  vector  d  there  exist  unique  vectors  eA, 
e,  and  dA  satisfying  the  loop  equations.  This  implies  that  the  inverse  of  I-  exists; 
otherwise,  there  is  either  no  solution  to  the  loop  equations  or  there  are  an  infinite  number 
of  solutions.  When  the  inverse  exists 


e  =  Z,,(M,A2)d 

Ij(M,A2);=  Mn  +  M12A2(I  -  M^AJ1  M21 
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Figure  2-6  Nominal  and  Robust  Performance  specifications  Diagram 
where  L/M.Aj)  is  called  a  lower  linear  fractional  transformation.  In  order  to  analyze  the 
performance  specifications,  a  fictitious  block  is  created  between  the  input  d  and  the  output 
e.  Figure  2-6  shows  the  new  structure.  The  set  of  all  allowable  blocks  is  defined  as 

B^^eA:  d(A,)£l} 

In  this  formulation  the  matrix  Mu=  L/M.O)  may  be  thought  of  as  the  nominal  map  and 
A 2e  B2  viewed  as  a  norm  bounded  perturbation  from  an  allowable  perturbation  class,  Aj. 
The  matrices  M,2,  Mg,,  and  and  the  formula  L(M,  )  reflect  prior  knowledge  on  how 
the  unknown  perturbation  (Ag)  affects  the  nominal  map,  Mn.  In  this  case  L(M,-)  is  related 
to  as  defined  earlier.  This  type  of  uncertainty,  called  linear  fractional,  is  natural 

for  many  control  problems,  and  encompasses  many  other  special  cases  considered  by 
researchers  in  robust  control  and  matrix  perturbation  theory.  The  constant  matrix  problem 
to  solve  is:  determine  whether  the  LFT  is  well  posed  for  all  A2  in  B2  and,  if  so,  determine 
how  "large"  L/M,Ag)  can  get  for  Age  B2.  Define  a  new  structure  A  as 

A:={o  A0jA'eWA’ 


2-10 


Now  there  are  three  structures  with  respect  to  which  (i  can  be  computed.  They  are  as 
follows:  p,(.)  is  with  respect  to  A,,  is  with  respect  to  Aj,  and  |iA(.)  is  with  respect 

to  A.  In  view  of  this,  |X1(Mn),  and  ^A  (M)  are  all  defined.  Theorem  4.2  from 

[PD93]  says:  The  linear  fractional  transformation  M A2)  is  well  posed  for  all  Ajg  B2  if 
and  only  if  1-  As  the  perturbation  Aj  deviates  from  zero,  the  matrix  L/M^) 

deviates  from  Mn.  The  range  of  values  that  p.jfL/M.Aj))  takes  on  is  intimately  related  to 
|aa(M),  as  follows: 


This  relationship  is  known  as  the  Main  Loop  Theorem  [PD93]. 


^2(^22)  ^  1 


A,6B, 


233  Structured  Singular  Value  Synthesis 

The  ^-synthesis  problem  is  described  by  the  attempt  to  find  a  controller  K(s)  that 
minimizes  an  upper  bound  on  the  structured  singular  value, 


inf  inf  supd[DM(K)D1] 

K  mhiHzmg  DeD  oo 

where  M(K)  is  the  closed  loop  transfer  function.  One  way  to  solve  this  problem  is  called 
D-K  iteration;  it  calls  for  alternately  minimizing  sup  C(DM(K)D'1 )  for  either  K  or  D 
while  holding  the  other  constant.  First  the  controller  synthesis  problem  is  solved  using 
design  on  the  nominal  design  model  (nominal  plant  plus  weighting  functions);  i.e,  Dsl. 
p-analysis  is  then  performed  on  the  closed  loop  transfer  function  M(K),  producing  values 
of  the  D  scaling  matrices  at  each  frequency.  The  resulting  frequency  response  data  is  fit 
with  an  invertable,  stable,  minimum  phase  transfer  function  which  becomes  part  of  the 
nominal  synthesis  structure.  With  D  fixed,  the  controller  synthesis  problem  is  again  solved 
by  performing  an  H„  design  on  the  augmented  system.  The  D-K  iterations  are  continued 
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until  a  satisfactory  controller  is  found  or  a  minimum  is  reached.  The  resulting  controller 
order  is  the  order  of  the  design  plant  and  weighting  matrices,  as  well  as  the  order  of  the  D- 
scale  transfer  function  fits  [ABSB92].  MatLab™  provides  a  p-toolbox  that  will  be  used  in 
this  thesis  which  performs  this  D-K  iteration. 

1.4  Nominal  Performance  .Robust  Stability  and  Robust 

Performance  tests 

This  section  presents  the  tests  for  Nominal  Performance  (NP),  Robust  Stability  (RS), 
and  Robust  Performance  (RP)  for  the  system  of  Figure  2-6.  Depending  on  the  type  of 
perturbation  (structured  or  unstructured),  the  infinity  norm  test  is  conservative  and 
^.-analysis  is  required  as  shown  in  Table  2-2. 


Table  2-2  Test  for  NP,  RS,  and  RP  Prom  [Doy85i 


Perturbation 

Stability  Test 

Performance  Test 

A=0 

No  C+  poles 

l,Mn>L^l  (NP) 

d(A)  <  1 

■MaA-si  (RS) 

P(M>S1  (RP) 

Ae  BA 

l*(Mn)£l  (RS) 

P(M><a  (RP) 

Table  2-2  summarizes  the  objectives  of  optimization  and  p-  synthesis.  Notice  that 

this  table  does  not  mention  any  test  using  the  two  norm,  and  the  objective  in  this  table  is 
only  to  minimize  II  .  Il„  or  p(.).  This  means  that  performance  based  on  white  Gaussian 
noise  inputs  is  not  accounted  for.  This  is  the  true  objective  of  the  mixed  H2/H„,  control 
design  problem:  to  address  the  tests  for  II .  Il„  and  provide  a  low  II .  as  will  be  seen  in 
the  next  chapters. 
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One  way  to  measure  the  robustness  of  a  system  is  to  calculate  the  Vector  Gain  Margin 
(VGM)  and  Vector  Phase  Margin  (VPM)  at  the  input  and  output  of  the  plant  (MIMO). 
The  VGM  and  VPM  tell  us  how  much  the  system  can  tolerate  a  change  in  gain  and  phase 
before  it  goes  unstable.  The  VGM  and  VPM  using  the  complementary  sensitivity  function 
is  measured  by 


r((0)  =  — : — 

o(TO'<o)) 

where  T  can  be  at  the  input  if  we  are  looking  at  the  input  margins  or  at  the  output  of  the 
plant  if  we  are  looking  at  the  output  margins.  Whichever  point  we  are  looking  at,  the 
general  formulas  are 

VGMT  =[1-^,1+^]  where  =  inf  r(co) 

<D€K 

and 

VFMt  =  [-0,+6]  where  0  =  2sm"^-^-j 
The  VGM  and  VPM  using  the  sensitivity  function  are  defined  through 


/>(«)= - - - 

d(S(yco)) 

where  S  can  be  at  the  input  if  we  are  looking  at  the  input  margins  or  at  the  output  of  the 
plant  if  we  are  looking  at  the  output  margins.  Whichever  point  we  are  looking  at,  the 
general  fi  rmulas  are 


VGMS  = 


_ 1 _ 1_J 

.l  +  Pmin’i-Pnm. 


where  pmjD  =  inf  r(co) 

ooeK 
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and 


VPMS  =  [-6,  +0]  where  8  =  2  an 

Since  VGMT,  VGMS,  VPMT ,  and  VPMS  are  all  important,  this  thesis  will  compute  all  of 
them  in  order  to  evaluate  the  level  of  robustness  at  the  input  and  output  of  the  plant.  For 
more  details  see  [Dai90]. 
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Constraint 


This  section  presents  the  mixed  optimization  developed  by  [Rid91].  Mixed 

Hg/HL  optimization  is  a  nonconsevative  tool  that  trades  between  Hj  and  H„  objectives. 
The  goal  of  die  mixed  problem  is  to  find  a  stabilizing  compensator  K(s)  that  achieves 

..  M.Jfarl,  .subject to  HT^IL^y 

where  Tw  and  T^  can  be  defined  to  be  independent  of  each  other.  Figure  3-1  shows  the 


Mode  diagram. 


e 

z 


Figure  3-1  Mixed  Design  Diagram 
The  state  space  matrices  are: 


X 

=  Ax 

+ 

Bdd 

+ 

Bww 

+ 

B0u 

e 

=  Cex 

+ 

Dgjd 

+ 

DewW 

+ 

z 

=  Cxx 

+ 

D^d 

+ 

+ 

y 

=  Cyx 

+ 

D^d 

+ 

DywW 

+ 
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3.1  Nopsingular  H™  Constraint 


The  following  assumptions  are  made  on  the  state  space  matrices: 
(i)  D^j  =  0  ;  (ii)  Dw  =  0  ;<m)D„«0 

(iv)  (A,  Bu)  stabilizable  and  ( Cy,  A)  detectable 

(v)  DT#nDeiI  full  rank;  DJldDT)d  full  rank 

(vi)  DT^D,,,  full  rank ;  DywDTyw  full  rank 


|"A  -  jtiJl  Bn  "I 

(vii)  has  full  column  rank  for  all  co 

L  Ci  DaJ 

rA-j©i  Bwi 

(viii)  I  I  has  full  row  rank  for  all  co 

|_  Cy  DywJ 

fA-jfi)I  Bn"! 

(ix)  ^  I  has  full  column  rank  for  all  to 

P A  -  j0£>  I  Bal 

(x)  I  _  has  full  row  rank  for  all  ca 

Cy  DydJ 

These  conditions  are  the  union  of  the  Bj  assumptions  and  the  1C  assumptions,  except  that 
the  controller,  K(s),  for  the  mixed  problem  must  be  strictly  proper  in  order  to  guarantee  a 
finite  two-norm  for  Tw.  The  state  space  matrices  for  K(s)  are: 


xc  =  Acxc  +  Bcy 
u  =  Ccxc  ;  Dc  =0 

and  the  closed-loop  matrices  are: 

x  =  Ax  +  Ba  d  +  Bw  w 

e  =  Cex  +  D^d  +  D"  w 

z  =  Ctx  +  D^d  +  Dww 
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=  0  and  =  0 


andTw  can  be  written  as 

Ted=Ce(s/-A)-,Bd 


The  following  definitions  will  be  made: 


irf  |Tj 

KMfaibans" 


T„  =C,(s/-A)'B„ 


K^.  s  the  unique  K(s)  that  makes  ,|TIwll2  =  a0 


Y  a  II  Trf  IL  when  K(s)  = 

*  a  K(s)  that  solves  the  mixed  Hj/H^  problem  for  some  y 
/  *  HTadILwhenK(s)  =  K1Ilix 

a*  =IITwrll2whenK(s)  =  Kmix 


Theorem  4.1.1  from  [Rid91]  says: 

Thcoran3,l :  Let  (Ac,  Bc,  CJ  be  given  and  assume  there  exists  a  QL  =  Ql  ^  0  satisfying 

AQ-  +  Q.A7  +  r'Q-C'.CJL  +  B„B]  =  0  (••) 

The  following  are  equivalent: 

i)  (A,Ba)  is  stabilizable 


ii)  A  is  stable 


Moreover,  if  i)  -  ii)  hold,  the  following  are  true: 


iii)  the  transfer  function  T^  satisfies 


iv)  the  two  norm  of  the  transfer  function  Tw  is  given  by 

|T„g  =  «CiaC'ITJ  =  tr[G2QTCJ 


where  Q2  =  Q2  >  0  is  the  solution  to  the  Lyapunov  equation 


AQ2+Q2AT  +  BwBl=0 
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v)  all  real  symmetric  solutions  to  (**)  are  positive  semidefinite 

vi)  there  exists  a  unique  minimal  solution  to  (**)  in  the  class  of  real  symmetric  solutions 
vn)fl.  is  the  minimal  solution  to  (**)  iff 

Rea,(A  +  Y'JQ-CjCt)]<0  V< 

A+y  -*Q.  ,CeTCe  is  stable,  where  <2L  is  the  minimal  solution  to 

(**) 

Proof:  See  Theorem  4.1.1,  [Rid91]. 

Using  Theorem  3.1,  the  mixed  problem  can  be  restated  as: 

Find  a  strictly  proper  controller  K(s)  that  minimizes  the  index 

J(At,Bc>Cc)  =  tr(&C,TC,) 

where  Q2  is  the  teal,  symmetric,  positive  semidefinite  solution  to 

AQ2+Q2AT  +  B„Bl=0 

and  such  that 

AQ. + <LAr  +  r2Q-CjC,a.  +  Bfil = 0 

has  a  real  symmetric  positive  semidefinite  solution.  To  solve  this  minimization  problem 
with  two  equality  constraints,  a  Lagrange  multiplier  approach  is  used: 

L = tr(Q2C,TCJ+  tr  ([AQ2  +  Q2Ar  +  BWB?]X) 

+tr  UQ.+a.AT  +Y“2e.CtTC,a.  +  BdBdT  in 
where  X  and  Y  are  the  Lagrange  multiplier  matrices.  The  necessary  conditions  for  a 
minimum  are  given  by  [Rid91].  Conclusions  from  these  conditions  are: 

(i)  Nc  mixed  solution  exists  for  7  <  y 

(ii)  The  mixed  solution  comes  from  seven  first  order  necessary  conditions,  which  are 
highly  coupled  and  nonlinear. 

(iii)  For  y<y<y,  neutrally  stabilizing  ARE  solutions  are  required,  and  y*  =  7. 
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(iv)  Far  y  2a  y ,  is  the  unique  mixed  solution. 


3.2  Singular  H„  C 


i]  i  Ki 


This  section  presents  the  mixed  optimization  problem  with  a  singular  H. 

constraint,  developed  by  [WR94a]. 


3.2.1  General  Formulation 

Here,  mured  H2/Hw  optimization  will  be  extended  to  handle  a  (possibly  non-strictly 
proper)  singular  H„  constraint.  Assume  the  plant  P  contains  the  H2  and  the  H„  designs. 
The  individual  Hj  and  designs  can  be  represented  as  two  independent  systems 


A2 

Bw 

b„; 

’A. 

Bd 

B»_‘ 

It 

ft- 

C, 

P-  = 

C« 

D. 

C,. 

Dy. 

D,» 

-Cr- 

V 

where 


design  ]  » 


c,-  =[c 


pint 


(3-1) 


The  objective  for  the  mixed  case  is  to  find  a  stabilizing  compensator  K(s)  that  achieves 

k.SL.I'1’’*!*  • subject  10  M- 5  Y  (3_2) 


where 

T«=C,(sI-Ay'Ba  +  D^  ;  T„  =Ct(sI-AT'Bw  +D„  (3-3) 

are  the  closed  loop  transfer  functions  from  w  to  z  and  d  to  e,  respectively.  The  following 
assumptions  are  made  in  die  state  space  matrices: 
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(i)  Dw  =  0  ;(ii)Dyu  =  0 

(iii)  ( A2 ,  BUj )  stabilizable  ami  (Cy,A2)  detectable 

(iv)  (A„,B„  )  stabilizable  and  (Cy  ,A„.)  detectable 

(v)  DT^D,,,  full  rank ;  DywDTyw  full  rank 


(vi) 


A2  -  j(Ol 


has  full  column  rank  for  all  ct> 


[A2-j(Ol  B1 

(vii)  I  r  n  has  full  row  rank  for  all  to 

Notice  that  D^,  is  not  restricted  to  zero  and  no  assumptions  are  made  as  to  the  ranks  of 
D",  and  Dyd ;  this  means  that  a  singular  H„  design  can  be  allowed.  For  the  mixed  problem, 

K(s)  must  be  strictly  proper  in  order  to  guarantee  a  finite  two-norm  for  T^.  The  state 
space  matrices  for  K(s)  are: 


x«  =  Acxc+Bcu 


u  =  Ccxc 


;D  =o 


(3-4) 


and  the  closed-loop  matrices  are: 


X2  =  A2x2  +  J3ww 
z  =  Czx2 


where 


x-=A„x„  +  Bd  d 

e  =  Ccx„  +  D^d 


(3-6) 


4.= 


,  B  C 

L  e  y_ 


Bu.Cc‘ 

Ac  . 


(3-7) 
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(3-8) 


Ct=[Cz  D.CJ;  Ce=[Ce  DmCc]  (3-9) 

Ad=D«i  (3-10) 

Using  the  definitions  from  Section  3.1,  the  mixed  H^/H^  problem  is  now  to  find  a 

controller  K(s)  such  that: 
i.  A2  and  A„  arc  stable 

ii- |TedL  -7 for  7^7 

iii.  II T,^  llj  is  minimized. 

Now  Theorem  3.1  can  be  restated  as  follows: 

Theorem  3.2:  Let  (Ac,  Bc,  CJ  be  given  and  assume  there  exists  a  solution 
(3L  =  <2l  >  0  satisfying 

K.Q~  +  Q-Al  +  (Q.CJ  +  Bfil  )K'(Q_Cl  +  Bfil  )T  +  BaB]  =0  (3-11) 

where  R  —  (y2I  -  )  >  0.  Then,  the  following  are  equivalent: 

i)  ( A. ,  )  is  stabilizable 

ii)  A„  is  stable. 

Moreover,  if  i)  -  ii)  hold,  the  following  are  true: 

a»rU.SY 

iv)  the  two  norm  of  the  transfer  function  Tw  is  given  by 

lT„|g=trlCie2C,I]  =  tr((22CjCJ 
where  Q2-Q2  -  0  is  the  solution  to  the  Lyapunov  equation 

AtQt+Qrf  +  BX^O 

v)  all  real  symmetric  solutions  to  (3-1 1)  are  positive  semidefinite 

vi)  there  exists  a  unique  minimal  solution  to  (3-1 1)  in  the  class  of  real  symmetric 
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solutions 

vii)  £?„  is  the  minimal  solution  to  (3-11)  iff 

Re[X.J(A_  +  BdD>,Ct  +  <2.CIT«IC,)]<0 

viii)  |T„|_  <  7  iff  (/4_  +  BdD^jR'C^  +  Q_CtTfi‘CJ  is  stable,  where  is  the 

minimal  solution  to  (3-11) 

Proof:  See  Theorem  3  [WR94a]. 

Using  Theorem  3.2,  the  mixed  case  can  be  restated  as: 

Find  a  strictly  proper  controller  K(s)  that  minimizes  the  index 

J(Ac,Bc,CJ  =  tr(ftC,TC,)  (342) 

where  Q2  is  the  real,  symmetric,  positive  semidefinite  solution  to 

AQ2+Q2AJ+BwBl=0  (3-13) 

and  such  that 

A.O.  +  G.AI  +  (e.C;  +  BJo;)/r1(e.C.T  +  BJDi)T  +  BX  =0  (344) 

(with  R>0)  has  a  real  symmetric  positive  semidefinite  solution.  The  Lagrangian  for  this 
problem  becomes 

L  =  ir[Q2CjC,] +tr{[A2Q2  +  QX  +  B„Bj  ]X} 

+"■([  4.0.  +  (LAl  +  (O.CJ  +  Bfil  )R-'  (QjCj  +  Bfil  )T  +  BX  in 

(3-15) 

where  X  and  Y  are  symmetric  Lagrange  multiplier  matrices.  The  resulting  first  order 
necessary  conditions  have  not  been  solved  analytically  but  do  provide  some  insight  into 
the  nature  of  the  solution.  In  particular,  the  condition 


= (A. + BXXc. + o.c.T/r'cjTr 
+  Y(A. + bXxK + Q.cX'cj= o 


(346) 
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implies  that  either  Y  =  0  or  (A„  +  BdD^R'lCe  +  QjCjR  xCc)  is  neutrally  stable.  The 
former  condition  means  the  solution  is  off  the  boundary  of  the  H„.  constraint,  and  the  latter 
solution  implies  the  solution  lies  on  the  boundary  of  the  constraint  and  £?»  is  the 
neutrally  stabilizing  solution  of  the  Riccati  equation.  The  necessary  conditions  for  a 
minimum  are  given  by  [WR94a].  However,  it  is  easy  to  show  that 

(i)  No  mixed  solution  exists  for  y  <  y 

(ii)  For  Y  <  y  <  y,  neutrally  stabilizing  ARE  solutions  are  required,  and  y*  =  y. 

(iii) Fory^y,Kmil  =  K2opt. 


UU  Numerical  Solution 

Walker  and  Ridgely  [WR94a]  developed  a  numerical  method  for  synthesizing  a  family 
of  general  mixed  Hj/H,,,  controllers  which  can  handle  a  proper,  singular  H„  constraint.  In 

the  single  constraint  mixed  problem,  for  y  <  y  <  y,  the  solution  to  the  mixed  problem 

must  lie  on  the  boundary  of  the  constraint.  Further,  for  y  <  y  <  y,  a*  is  a  monotonically 

decreasing  function  of  y.  Finally,  for  y  <  y ,  no  solution  exists.  This  is  shown  graphically 

in  Figure  3-2.  The  numerical  method  for  solving  the  mixed  problem  was  motivated  by 
Figure  3-2.  Since  the  optimal  H2  controller  is  relatively  easy  to  calculate  and  it  provides  a 
point  on  the  desired  curve,  it  was  selected  as  the  initial  controller.  The  optimal  H2/ 

curve  is  generated  by  starting  at  the  optimal  H2  controller  and  stepping  along  the  a  versus 
y  curve  by  reducing  y  from  y  to  y  by  increments. 
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Y  7=  /  Y 

II T  II 

ed  od 


Figure  3-2  Hj/H.  Boundary  plot 


Applying  the  results  from  the  previous  section,  it  is  seen  that  the  optimal  mixed  Rj/H.. 
controller  for  a  fixed  7  will  have  the  property  that  11  Ted  II-  =  Y-  This  suggested  a  penalty 


function  approach  to  the  problem.  Consider  die  following  performance  index 


/,=fr4+x(fU-Yr 


(3-17) 


where  A.  is  a  penalty  on  the  error  between  the  desired  7  and  the  infinity-norm  of  the 


transfer  function  T*,.  Define  the  vector  X  as 


*-[<  -  <  bl  -  bl  <  -  <]  <3-|8> 

where  ac  ,  bc  and  CCj  are  the  columns  of  Ac,  Bc  and  Cc,  respectively.  The  first  order 


necessary  conditions  for  to  be  a  minimum  are 


a 

=  0  for  i=l,...,(nxn+nxp+nxm) 
OX; 


*r.£  , 

dx,  dx; 


(3-19) 


(3-20) 
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where  Xj  are  the  elements  of  X  [WR94a].  The  first  term  on  the  right  hand  side  of  (3-20) 
can  be  solved  analytically  using  the  results  of  the  previous  section.  The  second  term, 
however,  represents  complex  matrix  relations  and  is  evaluated  analytically  using  the 
sensitivity  of  the  norm  developed  by  [GL93].  Assuming  the  maximum  singular  value 
of  has  a  single  peak  for  coe  then  the  derivative  of  the  infinity-norm  can  be  written 
as 


afc.L  _ 


dx. 


V, 

.  V  dx  Ax„n  m 


(3-21) 


where  ux  and  vt  are  the  singular  vectors  associated  with  T^,  c^,  is  the  frequency  where  the 
maximum  singular  value  reaches  its  peak  value,  is  the  nominal  X  vector,  and 
denotes  the  real  part  The  derivative  of  the  transfer  function  can  be  determined  from 

^ed  (®o  )| 


dx 


Xmm 


=  -7-(C.0'“-A.)1Bi  +  Dtf)j 

dX  \xnom 


(3-22) 


The  second  term  on  the  right  hand  side  of  equation  3-20  can  now  be  written  as 


(3-23) 


A  DFP  approach  similar  to  the  algorithm  described  in  [RMV92]  is  used  to  minimize  the 
performance  index.  The  basic  algorithm  is  as  follows: 

1.  Compute  the  optimal  Hj  controller  and  set  up  the  initial  X  vector 

2.  Compute  y  and  set  y  =  y 

3.  Decrement  y 

4.  Perform  DFP  search  over  X  vector  space  for  minimum  Jy 

5.  Store  resulting  controller  and  repeat  from  step  3. 


Initially,  the  algorithm  can  be  ran  with  loose  tolerances  on  the  DFP  search  to  define  the 
desired  a  versus  y  curve,  then  the  convergence  tolerances  can  be  tightened  and  a 
particular  point  can  be  refined  to  desired  accuracy.  In  addition,  this  algorithm  can  be 
applied  from  any  initial  condition,  not  just  the  optimal  Hj  controller,  by  substituting  the 
appropriate  initial  X  vector  and  y.  Finally,  the  X  vector  was  defined  with  a  fully  populated 
state  space  form;  by  using  canonical  forms,  the  number  of  variables  can  be  reduced. 
However,  there  are  drawbacks  to  using  some  canonical  forms  such  as  the  controllability 
canonical  form  due  to  numerical  instability.  The  modal  canonical  form  has  been  used 
successfully  to  reduce  the  parameter  space.  The  numerical  solution  for  the  mixed  problem 
with  a  singular  constraint  can  be  extended  to  allow  multiple  HM  constraints  as  will  be 

shown  in  the  next  chapter. 
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Hoo  Constraints 


Two  major  goals  in  a  control  design  are  to  design  controllers  which  yield  Nominal 
Performance  (NP)  and  Robust  Stability  (RS).  These  can  be  represented  as 

II  WjS  IL  <;  1  for  Nominal  Performance  (NP)  (tracking) 

II W2T IL  ^  1  for  Robust  Stability  (RS) 

(Multiplicative  perturbation  (1+  AW2)G;  II A  ll„<  1) 
Using  this  perturbation  model  and  the  NP  condition,  [DFT92]  defines  the  Robust 
Performance  (RP)  condition  for  a  SISO  system  as 


IIWjTIL^l  and 


-  WiS-  <i 

i  -|w2t|  . 


which  is  also  given  by 

II  IWiSI  +  IW2TIIU<1 

This  formulation  is  often  solved  with  a  mixed-sensitivity  approach,  which  penalizes  both 
Sensitivity  (S)  and  Complementary  Sensitivity  (T),  as 


This  mixed  sensitivity  cost  function  is  required  to  satisfy 

frJ.<i/V2 


in  order  to  have  RP.  In  the  author's  opinion,  this  method  is  conservative,  because  the 
designer  has  no  control  over  the  trade-off  between  RS  and  NP.  This  chapter  presents  a 
nonconservative  method.  Recall  the  conditions  for  NP  and  RS  at  the  start  of  this  chapter. 
Both  H_  conditions  are  solved  independently  using  mixed  Hj/H.  optimization  with 
multiple  H,  constraints.  For  a  MIMO  system,  the  objective  is  to  achieve  RS  to  certain 
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perturbations  and  NP  at  possibly  a  different  location  in  the  loop.  Therefore,  this  new 
technique  will  permit  the  exploration  of  these  different  objectives  independently. 

4.1  Development  of  Multiple  Constraints 

This  section  presents  the  mixed  H2/H„  optimization  problem  with  multiple  H„. 
constraints.  Mixed  H^/H.  optimization  is  a  nonconsevative  tool  that  trades  between  Hj 
optimization  and  multiple  constraints.  Consider  the  system  in  Figure  4-1,  where  dj, 
i=l,...,n„  are  of  bounded  energy  ( II  dj  ll2  £1)  and  w  is  of  bounded  spectrum.  The  transfer 

function  P  is  the  underlying  plant  with  all  weights  associated  with  the  problem  absorbed. 
It  is  assumed,  in  general,  that  there  is  no  relationship  between  eit  i=l,...,n„  and  z  or  dj  and 

w.  The  input  w  is  unit  intensity  zero-mean,  white -Gaussian  noise  and  the  inputs  dj  are  of 
bounded  energy.  In  general,  the  state  space  of  P  is  formed  by  wrapping  the  weights  from 
an  H2  problem  from  w  to  z  and  the  weights  of  the  H_  problems  from  dj  to  ej  around  the 

basic  system  resulting  in  the  augmented  plant 


Figure  4-1  Mixed  Hj/H..  with  Multiple  H„  Constraints  Design  Diagram 
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The  plant  P  contains  the  Hj  design  and  the  designs.  The  individual  H2  and 
problems  can  be  represented  as  different  systems,  where 


a2 

Bw 

Bu  ' 
“2 

A„ 

t 

< 

P2  = 

c. 

Dm 

zn 

P-  = 

D<A 

Den 

CjU 

C„ 

v 

K 

D* 

D„ 

i  l,...9nw  (4-2) 


where 


®ptam 
®Hwde«ign_ 


=[Cpa«  ;  Cy_,  =[Cpbint 

The  objective  for  the  mixed  problem  is  to  find  a  stabilizing  compensator  K(s)  that  achieves 

jsjru  •  ***** ®  M.  * i= 1 n-i  (4-3> 

where 


Twr=Cl(sI-A2)-,Bw+Dw 

T^C^sI-A^B^+D*, 

are  the  closed  loop  transfer  functions  from  w  to  z  and  dj  to  ej,  respectively. 


(4-4) 
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The  following  assumptions  are  made  in  the  state  space  matrices: 

(i)  Dw  =  0  ;(ii)Dp=0 

(iii)  (Aj.B^)  stabilizable  and  (Cyj,A2) detectable 

(iv)  (A^.B^)  stabilizable  and  (C^.A^)  detectable  for  alii 

(v)  DT  JD,,,  foil  rank ;  DywDTyw  full  rank 


(viii) 

(ix) 


A2  -joal 

.  cz 

A2  -  j(0l 


has  full  column  rank  for  all  co 


has  foil  row  rank  for  all  to 


Notice  that  the  Ded|  are  not  restricted  to  zero  and  no  assumptions  are  made  as  to  the 
ranks  of  De  u  and  .  This  means  that  singular  H_  constraints  can  be  allowed.  The 

controller  K(s),  for  the  mixed  problem,  must  be  strictly  proper  in  order  to  guarantee  a 
finite  two-norm  far  T„.  The  state  space  matrices  for  K(s)  are: 


Xc  =  Acx  +  Bcy 
u  =  Ccx  ;  Dc  =0 

and  die  closed-loop  matrices  are: 


X2  —  A2X2  +  Bw  w 
z  =Czx2 


x-i=Amixmi+Biidi 

ei  =  + 


(4-5) 


(4-6) 
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what 


°x1 

;  A*  = 

A- 

Ac  J 

(4-7) 

Bw  = 

‘  Bw  ' 
B  D 

_"c^yw_ 

=  V 

(4-8) 

S' 

II 

n 

H 

D  Cl;  C  - 

au  cj  »  e, 

[C,  Dc„Cc] 

(4-9) 

<P 

II 

O 

*. 

(4-10) 

The  following  new  definitions  are 

made: 

7,  s  |Xd,|_  when  K(s)  =  K2v 

a  solution  to  the  Hg/H.  problem  for  some  set  T 
Y>|T-,Lwho.K(s)=  K„ 
a3|T„{2wtenK(s)=  K* 

The  mixed  problem  is  now  to  find  a  controller  K(s)  such  that: 

i.  A2  and  are  stable  for  all  i 

ii.  jT^  j^  <  Yj  for  some  given  set  of  Yi  >  Y. 

iii.  II llg  is  minimized. 

Now  Theorem  3  2  can  be  extended  to  multiple  H_  constraints  as  follows: 

Theorem  4.1:  Let  (A^,  Bc,  CJ  be  given  and  assume  there  exists  a  solution 
CL,  =  Qli  -  0  satisfying 

AJL*  +  GLAi  +  (Q^Cl  +  BAiD\  )R;1((LCl  +  Bdp\  )T  +  =  0  (4-11) 
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for  all  i,  where  =  (Y f  /  -  )  >  0.  Then,  for  each  i  the  following  are  equivalent: 

0  is  stabilizahlc 

ii)  Am 4  is  stable. 

Moreover,  if  i)  -  ii)  hold,  the  following  are  true: 

m)  KL^iforalli 

iv)  the  two  norm  of  the  transfer  function  Tw  is  given  by 

lT„U=trICie2C,Tl  =  tr[e!C,TCJ 
where  Q2=Q2  ^  0  is  die  solution  to  the  Lyapunov  equation 

+  0 

v)  all  real  symmetric  solutions  to  (4-1 1)  are  positive  semidefinite  for  all  i 

vi)  there  exists  a  unique  minimal  solution  to  (4-1 1)  in  the  class  of  real  symmetric 
solutions  for  each  i 

vii)  j  are  the  minimal  solutions  to  (4-1 1)  iff 

+  B'DltfC',  +  CUC^'cJ]  £0  for  ally 
™i)|T-iL<Yi  iff  (Aj  +(L,ClR:'Cti)  is  stable,  where  g_,  are 

the  minimal  solutions  to  (4-11)  for  all  i 
Proof:  See  Theorem  1,  [UWR94]. 

Using  Theorem  4.1,  die  mixed  case  can  be  restated  as: 

Find  a  stricdy  proper  controller  K(s)  that  minimizes  the  index 

J(Ae,Bt,Ct)  =  tr(Q2CjCI)  (4-1? 

where  Q2  is  the  real,  symmetric,  positive  semidefinite  solution  to 

AQ2  +  Q2At  +  BwBl  =  0  (4-13 


and  such  that 


AJL* + OX  +  <QJZ + bad'  )K(QJo: + BaD\ )t + BX  =  0  (4-14) 


has  a  real  symmetric  positive  semidefinite  solution  for  all  i.  To  solve  this  minimization 
problem  with  many  equality  constraints,  a  Lagrange  multiplier  approach  is  used.  The 
Lagrangian  is 


L  =  ,r[Q2CjC,]  +  tr{lA!Q2+Q2Al  +  B„BTJX) 

+2^{[4_GLi  +  +  (O.C*  +  B,pl  )*,<ftllCj  +  B^ph  f  +  B,  Bl  J}') 


(4-15) 

where  X  and  Yi  are  symmetric  Lagrange  multiplier  matrices.  The  resulting  first  order 
necessary  conditions  have  not  been  solved  analytically  but  do  provide  some  insight  into 
the  nature  of  the  solution.  In  particular,  the  condition 


dL 

ao. 


=  (^  +  5d^^Cei  +<LiClR:'Ctl)TYi 
+*i<^ +BapZRl'Cti  +e,ctxc,)=o 


(4-16) 


implies  that  either  Yi  =0  or  (A^  +  +  QiJtC^R^Cz  )  is  neutrally  stable. 

The  former  condition  means  the  solution  is  off  the  boundary  of  the  corresponding 

constraint,  and  the  latter  solution  implies  the  solution  lies  on  the  boundary  of  the 
corresponding  H„  constraint  and  CLi  is  the  neutrally  stabilizing  solution  for  that  H„,  Riccati 

equation.  From  this,  it  is  not  hard  to  show  that: 

(i)  no  solution  to  the  mixed  problem  exists  if  Yi  <  y  for  any  i 

-i 

(ii)  the  solution  to  the  mixed  problem  is  the  Hj  optimal  compensates',  K^,  if  Yi  -  Yi 
fen- alii 

(iii)  if  neither  (i)  nor  (ii),  the  solution  to  the  mixed  problem  is  on  at  least  one  of  the  H. 
constraint  boundaries,  and  a  neutrally  stabilizing  ARE  solution  is  required. 
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Condition  (iii),  which  holds  for  any  "non-trivial"  mixed  problem,  poses  severe  numerical 
problems,  as  addressed  in  the  next  section. 

4J  Multiple  H~  Constraints:  Numerical  Methods 

Two  approaches  have  been  developed  to  compute  controllers  which  solve  the  mixed 
problem  for  multiple  constraints.  The  first  method,  called  the  Grid  Method, 
computes  the  set  of  controllers  which  satisfy  the  H„.  constraints  in  the  region  of  interest. 
This  is  accomplished  by  holding  all  but  one  constraint  constant  and  varying  the  remaining 
constraint.  The  second  method,  called  the  Direct  Method,  attempts  to  simultaneously 
reduce  all  H.  constraints.  For  the  remainder  of  the  discussion  it  will  be  assumed  that  there 
are  only  two  H.  constraints.  The  results  can  be  extended  as  necessary  to  handle  larger 
constraint  sets.  The  methods  are  based  on  the  performance  index 

A = Kt- +  a(M_ -y,f + Mb-l  -  yj  <«4> 

where  ^  are  penalties  on  the  error  between  die  desired  and  the  infinity-norm  of  the 
respective  transfer  function.  Note  that  this  requires  every  constraint  to  be  achieved 

with  equality,  which  is  not  necessarily  the  optimal  solution.  In  order  to  avoid  this,  the 
constraints  should  actually  be  treated  as  inequality  constraints,  which  requires  a 
constrained  optimization  method.  This  has  been  accomplished  using  Sequential  Quadratic 
Programming;  see  [Wal94].  In  this  thesis,  the  constraints  will  be  treated  as  equality 
constraints,  however.  Since  a  large  portion  of  the  "active”  region  will  be  mapped  out,  this 
poses  oily  a  small  restriction.  Furthermore,  as  it  has  been  shown  that  the  optimal  order 
problem  has  all  constraints  satisfied  at  equality  [WR94c],  the  controllers  found  here  are 

the  closest  fixed  order  controllers  in  a  two-norm  sense  to  the  optimal  (free  order) 

controllers.  The  resulting  numerical  optimization  is  basically  that  of  Section  3.2.2,  except 
with  additional  similar  terms. 
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The  grid  method  consists  of  solving  a  series  of  mixed  problems  by  holding  one  HM 
constraint  constant  and  reducing  the  second.  Once  the  optimal  curve  has  been 
determined,  the  first  constraint  is  decremented  and  the  process  is  repeated.  The  initial 
conditions  for  the  method  are  determined  by  solving  the  two  single  constraint  mixed 
Ha/HM  problems  to  define  the  region  of  interest.  The  process  results  in  a  grid  defined  by  a 

versus  y,  versus  y2.  The  resulting  grid  is  shown  in  Figure  4-2. 


Figure  4-2  Grid  Method 


The  introduction  to  this  chapter  suggested  this  method.  Since  the  design  objectives  are 
limited  to  a  specific  region,  one  approach  to  synthesizing  a  controller  would  be  to  reduce 
both  constraints  to  the  desired  level  without  computing  the  entire  grid  mentioned  in  the 

previous  approach.  Hie  direct  method  does  this  by  concurrently  reducing  the  constraints. 
The  process  used  in  this  approach  is  to  begin  at  the  optimal  H2  controller  and 
simultaneously  reduce  y,  and  y2  until  the  controller  is  found  which  meets  both  objectives. 
This  results  in  a  controller  of  fixed  order  which  meets  both  the  H„  constraints  and  has  the 

smallest  two-norm  for  the  H2  transfer  function.  Figure  4-3  shows  this  method.  Notice 
that  Oy  proper  selection  of  the  step  size  of  the  H„  constraints,  the  designer  can  select  a 

desired  direction.  Also  note  the  "hooks”  at  the  end  of  each  curve.  These  are  the  result  of 


Figure  4-3  Direct  Method 


Figure  4-4  Direct  Method  3D  curve  and  2D  projections 

the  trade  off  between  the  H_  constraints  encountered  as  approaches  y  .  Figure  4-3 

-i 

shows  that  the  resulting  curve  is  a  3D  curve  (JJT^  ||2  vs.  |T(d(  ||  vs.  ||Tedj  j|  ).  Therefore, 
for  any  3D  curve  there  are  three  projections.  These  2D  curves  are  the  j|2  vs.  ||Tcdi  j|^ 

curve,  the  IT^,  ||2  vs.  curve,  and  the  It^J  vs.  |jTedi  ||  curve.  This  is  shown  in 

Figure  4-4,  and  the  curves  are  denoted  as  curve  a,  b  and  c,  respectively. 

^.-Feasible  Solutions 

Assume  for  this  section  that  numerical  problems  in  computing  a  solution  do  not  exist 
There  are  boundaries  in  the  mixed  problem  where  feasible  solutions  do  not  exist.  These 
boundaries  are: 

(i)  No  controller  results  in  a<a0 

(ii)  No  mixed  solution  exists  for  y{  <  y  for  either  i=l  or  2 
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(iii)  For  certain  values  of  %  no  mixed  solution  exists  which  simultaneously  satisfies  all 
<x>-norms  constraints,  even  though  y{  >  y  ,  Vi 

-  i 

First,  define  three  planes: 

Plane  OL0:  the  plane  defined  by  letting  a=  a0  for  all  y,  i=lJ2.  All  a's  above  this  plane 
represent  a  subopdmal  solution  to  the  Hj  problem,  and  solutions  below  this  plane  are  not 
feasible. 

Plane  y  :  the  plane  defined  by  letting  y,=  y  ,  for  all  a,  y2-  All  yj  above  this  plane  are 
subqptimal  solutions  to  the  corresponding  H^j  design,  and  solutions  below  this  plane  are 
not  feasible. 

Plane  y  :  direct  analogy  of  plane  y  . 

-2  - 1 

Figure  4-5  shows  these  planes. 


Figure  4-5  HJK-i  feasible  planes 


Mixed  Hj/H^  design  is  a  tool  that  trades  among  the  Hj  design  and  the  different  H„ 
designs.  The  trade-offs  taken  two  at  a  time  are  examined  next. 


2  m  m  i  ■  tii  1 1  i  gf  i  m  h  i  n  i  m  j  1 1 4i  \  gg  m  n  r 


[Rid91]  showed  that  a*  is  a  monotonically  decreasing  function  of  ?  for  a  single 
constraint  mixed  problem.  Therefore,  a*  is  now  a  monotonically  decreasing  function  of  y,. 
If  the  optimal  a  versus  curve  is  computed,  then  unfeasible  solutions  lie  below  this  curve, 
and  any  solution  above  this  curve  is  feasible  but  suboptimal.  Graphically,  this  is  shown  in 
Figure  4-6.  The  numerical  method  computes  a  suboptimal  curve  that  is  close  to  the 
optimal  curve.  This  is  due  to  the  requirement  of  rinding  a  mixed  solution  numerically. 


Figure  4-6  Trade-off  among  Hj/H.,  (feasible  solutions) 
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First  consider  that  the  mixed  problem  has  only  |Tedi  as  a  constraint  (a  single 
constraint  problem).  In  this  case  the  mixed  problem  is  just  a  mixed  Hj/H..,  design.  For 
each  point  on  curve,  can  be  computed.  Next,  consider  that  the  mixed 

problem  has  only  |Todi  jj  as  a  constraint  (a  single  constraint  problem).  In  this  case,  the 

mixed  problem  is  just  a  mixed  Hj/H^  design.  For  each  point  on  Hj/H^  curve,  |Tedj  ||^ 

can  be  computed.  The  resulting  curves  from  the  two  different  designs  are  shown  in  Figure 
4-7.  Define,  the  Hj/H,.,  curve  as  the  optimal  curve  for  mixed  design,  and  the 

curve  as  the  optimal  curve  for  the  mixed  Hj/H^  design.  These  two  curves  define 
the  boundary  between  the  region  of  sub-optimal  solutions  and  the  region  of  "optimal" 
solutions  as  shown  in  Figure  4-7.  From  a  control  point  of  view,  we  are  interested  in  the 


Figure  4-7  Design  region  in  the  ||T(di  jj  /  ||Xedj  jj  plane 
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region  on  or  below  the  optimal  mixed  Hj/H^  design  curves;  that  is,  the  region  of  optimal 
solutions.  The  region  above  these  curves  is  not  of  interest,  as  shown  in  [WR94c].  Here, 
the  optimal  solution  "snaps  back"  to  the  optimal  single  constraint  curve,  and  thus  is 
subop  timal. 

Consider  now  that  the  mixed  problem  has  both  designs  as  constraints.  In  this  case, 
die  mixed  problem  is  a  mixed  H2/Ha.1/HM2  design.  There  exists  a  boundary  close  to  the  7 

—  i 

values  where,  below  this  boundary,  no  feasible  solutions  exist.  This  is  shown  graphically 
in  Figure  4-7.  This  boundary  is  difficult  to  find  analytically  or  numerically.  This  was  the 
region  iii)  alluded  to  at  the  beginning  of  section  4.3. 

Joining  the  planes  and  boundaries,  a  region  of  feasible  solutions  can  be  drawn,  as 
shown  in  Figure  4-8.  Regions  of  subop  timal  solutions  are  also  plotted.  Figure  4-9  shows 
a  surface  for  the  mixed  H2/H-,/Hw2  design.  On  this  surface  we  are  interested  in  the  solid 

checkered  region  that  corresponds  to  the  optimal  region,  especially  at  the  "knee"  where 
the  fs  are  close  to  the  optimal  values.  Sub-optimal  regions  (shaded  checkered)  are  also 
shown  in  Figure  4-9.  These  are  not  optimal  mixed  solutions,  since  their  values  of  7  are 
greater  than  those  for  the  optimal  curves  corresponding  to  the  mixed  F^/FLi  design  and 
the  mixed  H2/Ha>2  design,  respectively.  However,  these  subop  timal  regions  help  to 
visually  clarify  the  problem. 

The  next  chapter  will  present  a  SISO  example  as  an  introduction  to  this  new  synthesis 
method.  It  will  show  the  boundaries  that  were  discussed  in  this  chapter. 
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Figure  4-9  Hj/H„  Surface  of  Optimal  and  Suboptimal  solutions 
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This  chapter  illustrates  mixed  Hj/H*.  design  with  single  and  multiple  H_  constraints. 

The  numerical  method  is  that  developed  by  [WR94a],  which  permits  generalization  of  the 
H_  constraint,  (i.e,  and  /  or  DyjDj,  not  required  to  be  full  rank  and  #  0 

allowed).  Fen*  this  SISO  example,  the  objective  is  to  show  some  of  the  boundaries  and 
methods  discussed  in  the  last  two  chapters.  An  acceleration  command  following  design 
for  the  F-16  is  desired.  The  F-16  plant  consists  of  a  short  period  approximation  (a,  q),  a 
time  delay  (5)  [first  order  Pad6  approximation],  and  a  first  order  actuator  servo.  The  state 
space  matrices  are: 

'  -20  0  0  O' 

-0.188  -1.491  0.996  0 

A*=  -19.04  9.753  -0.096  0  (M) 

-4.367  35.264  -0.334  -40 

'20' 

0 

B,=  ;  C,  =[4.367  -35.264  0.334  80] ;  Dg  =[0]  (5-2) 

0m 

The  poles  and  zeros  of  the  plant  are: 

poles  X=[  -40.0000;  4.3535;  1.9025;  -20.0000] 
zeros  =  [  -1.2564  +1 1.9340i;  -1.2564  -1 1.9340i;  40.0000] 


Three  designs  are  produced:  an  Hj  desigr  a  mixed  Hj/H^  design  with  a  single  constraint, 
and  a  mixed  Hj/H^  design  with  multiple  constraints. 

5.2  H?  Design 

The  Hj  design  is  set  up  as  a  basic  LQG  problem,  as  shown  in  Figure  5-1.  In  the  H2 
design,  two  exogenous  inputs  w,  and  w2  enter  the  plant;  they  are  zero-mean  white 
Gaussian  noise  with  unit  intensity.  It  is  desired  to  minimize  the  energy  or  two  norm  of  the 
controlled  outputs  z,  and  The  weights  related  with  w,  ,w2 ,  Zj  and  Z2  are: 

Wind  disturbance  weight:  The  wind  disturbance  constitutes  an  exogenous  input  (w,).  It 
passes  through  a  coloring  filter  Wd  and  a  distribution  matrix  ¥,  where 

"  0  ' 

-1.491 
9.753 
_35.265 

Measurement  noise:  The  measurement  noise  is  represented  by  an  exogenous  input  (w2). 
w2  is  added  to  the  feedback  signal.  The  weight  for  w2is: 

W„  =  0.025 

Control  Usage:  The  weight  for  control  usage  z,  is  the  scalar 

Wc=1.0 

Normal  acceleration  output  weight:  It  represents  the  weight  on  the  normal  acceleration 
(NJ;  this  weight  is  chosen  to  be  the  scalar 

Wz=1.0 


wjW-*2!5L;  <F  = 

S  +  6.7 
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Therefore,  the  system  P2  is 


Figure  5-1 H2  Block  diagram  (F-16) 


_  rr. 

The  basic  conditions  that  are  checked  here  include  Dw  =  0,  Dyu  =  0,  full  rank, 

and  D*  Dffl  full  rank;  these  are  met  by  the  design  with  a  non-zero  Wn  and  Wc.  Therefore, 
the  design  diagram  is  properly  set  up.  Table  5-1  shows  the  results. 

_ Table  5-1 H,  Results 


a„  =  0.2530 

■mss 

s 

[-5.8609, 28.76861 

±57.6035 

T 

[-6.4183,  3.6504] 

±30.2811 

Although  the  objective  was  to  design  a  pure  regulator,  from  Table  5-1  we  see  that  the  H2 
controller  provides  good  margins.  The  VGM  and  VPM  are  based  on  the  magnitude  plots 
of  sensitivity  and  complementary  sensitivity  as  explained  in  Chapter  2,  Section  2.5.  Figure 
5-2  shows  the  magnitude  plot  of  sensitivity  and  complementary  sensitivity. 


Figure  5-2  Magnitude  of  Sensitivity  and  Complementary  Sensitivity  (dB) 
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The  magnitude  of  the  sensitivity  shows  that  the  H2  controller  attenuates  the  wind 
disturbance.  The  magnitude  of  complementary  sensitivity  represents  the  measurement 
noise  feedthrough  to  the  plant  output,  the  inverse  of  the  allowable  multiplicative 
uncertainty,  and  the  closed-loop  tracking  transfer  function.  Evidently  this  design  does  not 
provide  good  tracking  since  the  gain  is  above  0  dB  at  low  frequencies  and  rolls  off  too 
early. 


Two  mixed  designs  with  single  constraints  are  solved.  The  first  design 
(FtyHL,,)  represents  a  sensitivity  constraint  (T^  =W,S),  and  the  second 
represents  a  complementary  sensitivity  constraint  (Tedj  =W,T).  Two  objectives  are  set  up: 
the  first  objective  is  to  compute  the  infinity  norm  of  with  the  controller  obtained 
from  the  mixed  design,  and  the  second  objective  is  to  compute  the  infinity  norm 

of  Ted(  with  the  controller  obtained  from  the  mixed  Hj/H^  design. 


■m*  gvfTif!  1 !  B  ( fminTTfHTt  In  1  iT"  fui  ■>  j  WTWIi 


The  block  diagram  for  the  sensitivity  constraint  design  is  shown  in  Figure  5-3. 


Figure  5-3  Mixed  Block  Diagram  (Sensitivity  Constraint) 


5-5 


The  transfer  Auction  between  the  exogenous  output  et  and  the  exogenous  input  d|  is 
T-,  ,  and  is  given  by 

Tedi  =  WiS  (5-5) 


The  weight  for  sensitivity  is  a  low  pass  Alter  W,,  given  by 

100 


W,(s)« 


s  +  0.1 


The  objective  for  the  mixed  design  is 

subjecttolT^.L  <y, 

The  system  is 


(5-6) 


1 

> 

4, 

Ba., 

p 

II 

C, 

F*. 

D„ 

with  the  state  space  matrices  given  by 


(5-7) 


Ag  0  B| 

*1  l”Bg 

;  B,  » 

;  Bn  = 

_0  A.J  “•  Lb. 

■j  n-  Lb.. 

C„=[0  C,];  =[Di];  D=1„=[D.]  (5-8) 

C,_=[C,  0];  D^=[0];  D„=[0] 


Since  W,  is  a  strictly  proper  transfer  function,  D,  =  0,  and  therefore  and 

are  not  foil  rank  matrices.  Thus,  we  have  a  mixed  optimization  problem 
with  a  singular  constraint  The  performance  index  for  the  numerical  solution  is 

(5-9) 
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Figure  5-4  ^  vs.  curve 

Starting  from  the  optimal  H*  controller  and  stewing  along  the  a  vs.  7,  curve  by  reducing 
from  Yj  to  y  by  increments  produces  Figure  5-4. 


LOiiTu  MTlWTI 


iTmi  <  jp Tm 


Te(Jj  is  the  transfer  function  between  the  exogenous  output  e-j  and  the  exogenous  input 


d2  as  shown  in  Figure  5-5.  Here 


T*  =  WtT 


(5-10) 


with  the  weight  for  complementary  sensitivity  denoted  by  W,,  and  given  by 

1000  *  (s  +  0.01) 


(s  +  1000) 


Figure  5-5  Mixed  Hj/Hk2  Block  Diagram  (Complementary  Sensitivity  Constraint) 


The  objective  for  the  H2/H,„  2  design  is 


The  system  P„,2  is 


subjec"° 


(5-11) 


1 

> 

l 

B<, 

>C 

II 

c. 

^eadj 

D-,» 

c,., 

Dvd 

y“j 

where  the  state  space  matrices  are  given  by 


rAg 

°1 

;  b,=H 

B  =  [B‘l 

Lo 

A.J 

■  Loj 

-  Lb'J 

[o 

C.]; 

D«,d,  =[oj; 

DC,„=[D,] 

[C« 

o]; 

D*,=[0]; 

D„=[0] 

(5-12) 


(5-13) 
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Figure  5-6  fTJk  vs.  curve 

Note  that  is  zero,  which  implies  we  have  a  mixed  optimization  problem 

with  a  singular  constraint.  The  performance  index  for  the  numerical  solution  is 

Jt,=llT»C+^(|T«,L-rJ  (5-m) 

Starting  from  the  optimal  Hj  controller  and  stepping  along  the  a  vs.  y2  curve  by  reducing 
from  y2  to  Y2  by  increments  produces  Figure  5-6.  The  "ripples"  are  due  to  numerical 
inaccuracies,  and  due  to  the  monotonic  property,  could  be  "smoothed". 

5«3«3  Boundaries  on  the  mixed  l MIL*  aicfii cfc 

Using  the  controllers  K(s)  from  the  Hj/H^  design,  the  infinity  norm  of  T^  can  be 

computed,  which  represents  the  Complementary  Sensitivity  constraint.  Figure  5-7  shows 
the  trade  off  between  (weighted  Complementary  Sensitivity  constraint)  and  Tedi 

(weighted  Sensitivity  constraint).  Figure  5-7  shows  that  although  |Tedi  |j  is  being 
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Figure  5-7  IjT^J  vs.  curve  (Hj/H..,  design ) 

minimized,  there  isn't  a  big  change  in  |Tolj  J|  until  |Tedi  ||^  reaches  the  knee  in  the  a  vs.  7! 

curve.  Now  with  the  controllers  K(s)  from  the  H^HL2  design,  the  infinity  norm  of  T*.  is 

computed,  which  represents  the  Sensitivity  constraint.  Figure  5-8  shows  the  trade  off 
between  (weighted  Complementary  Sensitivity  constraint)  and  Tedi  (weighted 

Sensitivity  constraint).  Figure  5-8  shows  a  slightly  different  type  of  behavior  than  Figure 

5-7;  in  Figure  5-8  the  trade-off  in  KL  occurs  almost  immediately.  A  graphical 

interpretation  showing  the  relationship  between  the  two  norm  and  the  infinity  norms  is 
shown  in  Figure  5-9.  The  two  3D  curves  (solid  lines)  are  the  curve  for  the  H2/H„,j 
design  and  the  curve  for  the  Hj/H^  design.  The  dotted  lines  represent  the  projection  of 

the  3D  curve  for  the  design  on  the  vs.  ItJ^  vs.  |Todi  ||_ ,  and 

\\TJI  vs.  JtJ  planes  and  the  dashed  lines  represent  the  projection  of  the  3D  curve 
for  the  design  on  the  same  set  of  planes. 
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gure  5-8  **•  ML  curve  (H2/H..2  design ) 


gure  5-9  3D  lITJk  vs.  jTj^  vs.  curve  for 


H2/H„,  design  and  design 
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These  two  3D  curves  represent  the  boundaries  of  the  mixed  and  Hj/H „2  design  as 

mentioned  in  Chapter  4,  Section  4-3.  Unfortunately,  the  mixed  H2/H„  with  a  single  H,. 
constraint  problem  does  not  permit  both  H„  constraints  to  be  made  small  at  the  same  time, 

unless  they  are  wrapped  in  one  transfer  function.  The  next  section  presents  the  multiple 
constraint  results. 


5.4  Mixed  HVfiL.  Optimization  with_MultipteJL»  Constraints 

Now,  the  mixed  problem  is  a  mixed  H2/H0.1/H^2  design,  and  the  objective  is  to  find  a 
stabilizing  controller  K(s)  that  achieves 


inf  IItII  subject  to 

K  stabilizing 11  w|12 


(5-15) 


Again,  recall  that  both  inequality  constraints  will  be  treated  as  equality  constraints.  Thus, 
the  performance  index  for  the  numerical  method  is 

-hUC+MKL -r>)2 +^(IT-I  -TaJ*  f5-16* 

Two  approaches  are  used:  the  first  one  will  be  the  grid  method  and  the  second  one  the 
direct  method,  as  mentioned  in  Chapter  4. 


5,4»l£rrid  Method. 

The  results  of  applying  the  direct  method  (explained  in  Chapter  4,  Section  4.2.1)  are 
shown  in  Figures  5-10  and  5-11.  Figure  5-10  shows  the  |Xed2 1  vs.  ||Ted|  |  curves.  The 
lower  dotted  curve  represents  the  boundaries  for  both  H*,  constraints.  Also,  notice  that 

when  smaller  values  of  infinity  norms  are  reached ,  the  trade-off  between  the  infinity  norm 
of  the  constraints  starts  to  have  the  effect  discussed  in  Chapter  4,  Section  4.3.  Figure  5-11 
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Two_norm  Tzw 


shows  the  3D  surface.  This  surface  has  an  almost  flat  bottom;  therefore,  the  increase 
die  two  norm  starts  when  the  knee  of  the  individual  curves  is  reached. 


3-4-2  Direct  method 

The  direct  method,  explained  in  Chapter  4,  Section  4.2.2,  minimizes  both  H„ 

constraints  at  the  same  time.  By  proper  selection  of  the  steps  Ay,  and  Ay2  >  where 
Ayi  =  Yj  -  Yi  (for  i=  1,...  ,11*0  and  j=l,...,  number  of  steps),  the  designer  can  guide 

the  direction  of  minimization  to  the  desired  infinity  norms  ||T’odi  |  and  .  Figure 

5-12  shows  IT^IL  vs.  J|Tedi  |  for  four  different  directions.  The  starting  controller  was 


for  three  of  them  and  a  (taken  from  the  grid  method)  for  the  last  one. 


ITJ 
*1  * 

Figure  5-12  jjT^J  vs.  jr^l  Direct  Method 
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Figure  5-13  3D  plot  Direct  Method 
irrjb  vs.  M. 


Notice  that  this  numerical  method  permits  the  user  to  start  at  any  point;  this  means  the 
routine  can  start  with  any  stabilizing  controller.  Figure  5-13  shows  the  different  curves  in 
3D.  Finally,  Figure  5-14  and  5-15  show  the  final  results  of  both  the  grid  and  the  direct 
methods  combined. 
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S3  Conclusions  for  the  Multiple  H_.  Constraints  SISO  case 

This  new  optimization  technique  permits  minimization  of  the  two  norm  of  one  transfer 
function  subject  to  multiple  H_  constraints.  It  controls  the  trade  off  between  the  different 
constraints  and  the  H2  performance. 

Suppose  a  SISO  plant  has  a  multiplicative  perturbation  and  a  performance  requirement 
cm  sensitivity;  a  controller  K(s)  provides  Robust  Performance  for  that  plant  if  and  only  if 

II  I  WtSI  + 1 W2TI IL  <  1  (SISO) 

In  the  mixed  HJ  optimization  problem  with  multiple  H„  constraints  ,  the  two  H„ 
constraints  are  adjusted;  therefore,  it  does  not  have  frequency  specific  information. 
However,  if  we  define  the  following  requirements: 

II  W,S  ll„=  ft  <  1  and  occurs  at  any  CO  =  0^  (Nominal  Performance) 

and 

II W2T  IL=  ft  <1  and  occurs  at  any  co  =  cOf  (Robust  Stability) 

where 

y  <  ft  and  y  <  ft. 

“5  1 

then  two  cases  for  Robust  Performance  could  exist: 

Case  1.-  Hie  Robust  Performance  test  is  passed  iff 

ft  +  fc<l 

This  means  that  the  worst  case  for  a  Robust  Performance  test  is  that  II  W,S  ll—  and 
II W2T  IL  occur  at  the  same  frequency  (cog  =  ooj.). 

Case  2.-  If 

ft  +  ft  >  1 

the  Robust  Performance  test  must  be  applied  This  means  that  it  is  necessary  to  check 
frequency  information. 
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These  two  cases  relax  the  requirement  for  the  mixed  sensitivity  cost  function 
<11  T«d  II-  <  1/V2),  where 


because  now  the  designer  has  control  over  the  different  infinity  norms.  A  graphical 
interpretation  is  shown  in  Figure  5- 16. 


Nominal  Perform*™* 


Figure  5-16  Application  of  the  mixed  problem  with  multiple  H.  constraints 
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VI.  THE  HIMAT  PROBLEM:  A  M1MO  EXAMPLE 


For  a  MIMO  example,  the  HIMAT  problem  from  the  p-Tools  Manual  [Mat]  was 
selected.  The  HIMAT  vehicle  is  a  scaled-down,  remotely  piloted  vehicle  (RPV).  The 
design  example  will  consider  longitudinal  dynamics  only.  For  more  information  about  this 
problem  refer  to  the  p-Tools  Manual  [BDGPS91]  and  [SLH81].  The  HIMAT  problem 
will  be  designed  for  single  and  multiple  ^constraints. 


61  Problem  Stet  Up 

The  states  variables  of  the  plant  (HIMAT)  are: 
v  -  forward  speed 

a  -  angle  of  attack  (not  to  be  confused  with  a  from  the  Hj  design) 
q  -  rate  of  pitch 
0  -  pitch  angle 
Control  inputs  (u): 

5e  -  elevon  command 
Sc  -  canard  command 
The  variables  to  be  measured: 
a  and  6 
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The  state  space  matrices  are: 


'-0.0226  -36.6 

-18.90 

-32.1' 

0 

0 

0  -1.90 

0.983 

0 

Bg  = 

-0.414 

0 

0.0123  -11.7 

-2.63 

0 

-77.80 

22.40 

0  0 

1.0 

0 

0 

0 

JO  57.3000  0  0  1  Jo  ol 

[0  0  0  57. 3000 J  |_°  °J 

Short  period  roots  =  -2.2321  ±  3.3779i ;  Phugoid  roots  =  -0.0442  ±  0.2093i 

6.2  H*>  Design  requirement 

The  objective  of  the  Hj  design  is  to  develop  a  regulator  that  limits  the  white  noise 
feedthrough  to  the  angle  of  attack  and  pitch  angle  plant  outputs  (yg)  and  the  control 
usage  (u).  Figure  6- 1  shows  the  Bj  regulator  design  plant  with  the  weights  Wc  and  W2  on 
the  control  usage  and  states,  respectively.  Energy  from  the  white  noise  inputs,  w,  and  w2, 
will  be  minimized  with  respect  to  the  chosen  outputs,  z5  and  Zj,  by  the  compensator 
design. 

6.2.1  Weight  Selection 

Wind  disturbance  weight:  The  wind  disturbance  constitutes  an  exogenous  input  (w,). 
It  passes  through  'P  as  an  angle  of  attack  perturbation 

'F  =  [-36.6  -1.90  -1L70  of 


Figure  6-1  H2  Regulator  Diagram  (HIMAT) 

Measurement  i.uise:  The  measurement  noise  is  represented  by  an  exogenous  input  (Wj). 
w2  is  added  to  the  feedback  signal.  The  weight  for  w2  is: 

Wn=  0.1*Ip) 

Control  Usage:  The  weight  for  control  usage  Wc  is  chosen  as  the  identity  matrix 

Wc  =  I(2) 

Weighted  output:  This  is  the  weighted  angle  of  attack  and  pitch  angle;  here,  this  weight  is 
chosen  to  be  the  identity  matrix 

w,  =  iw 

The  state  space  for  the  Hj  design  plant  is 


A2 

Bw 

c 

D 

D 

z 

ZW 

2U 

D- 

(6-1) 
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and  the  state  space  matrices  are 


M  =  [  A,  lx,]  +  [«F  0][3']  +  [  Bg  Ju] 

ft]  ■  W-j  *  C  a:;]  *  &  « 

[y]  =  [  c.  W  *  [o Wn ♦  [oI„] 

Notice  that  the  "  0  "  represents  a  2ero  matrix  with  the  corresponding  dimensions  in  the 
different  state  space  matrices.  The  basic  conditions  that  are  checked  here  include  Dw  =  0, 
Dyu  =  0,  and  D^D^,  full  rank;  these  are  met  by  the  design  with  a  non-zero  Wn 

and  Wc. 


6*2,2.  HzJtesulte 


Table  6-1  shows  the  results. 


Table  6-1 H2 


Results  (H 


VGM  input  (dB) 

VPMi  (deg) 

VGM  output  (dB) 

VPMo  (deg) 

T 

[-0.1198,  0.1181] 

±0.7847 

[-11.7255,  4.8147] 

±43.4773 

S 

[-0.1181,  0.1198] 

±0.7847 

[-5.1530,  14.4199] 

±47.7754 

Although  the  objective  was  only  to  design  a  pure  regulator,  from  Table  6-1  we  see  that  the 
Hj  controller  provides  robustness  at  the  output  of  the  plant,  but  poor  margins  at  the  input 


of  the  plant.  The  VGM  and  VPM  are  based  on  the  magnitude  plots  of  sensitivity  and 
complementary  sensitivity.  Figure  6-2  shows  the  maximum  singular  values  of  the 
complementary  sensitivity  at  the  input  (Tj )  and  output  (TJ  of  the  plant.  This  represents 
the  measurement  noise  feedthrough  to  the  plant  input/output  and  the  inverse  of  the 
allowable  multiplicative  uncertainty  at  the  input/output  of  the  plant. 


Figure  6-2  Maximum  Singular  Values  of  Complementary  Sensitivity 


Figure  6-3  Maximum  Singular  Values  of  Sensitivity 
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Therefore,  provides  a  good  level  of  robustness  at  the  output  of  the  plant,  but  the 
system  is  susceptible  to  a  multiplicative  uncertainty  at  the  input  of  the  plant.  Figure  6-3 
shows  the  singular  values  of  input  sensitivity  (S;)  and  output  sensitivity  (SJ.  Notice  from 
Figure  6-3  that  the  input  sensitivity  is  minimized  for  good  wind  disturbance  rejection,  as 
seen  in  the  low  gain  in  S;.  The  problem  with  this  design  is  poor  performance  at  the  output 
of  the  plant,  as  seen  by  S0.  therefore  produces  a  system  that  is  weak  to  a 
multiplicative  perturbation  at  the  input  of  the  plant  and  has  bad  tracking  properties,  as 
shown  in  Figure  64. 

The  designs  to  be  addressed  are:  recover  the  margins  at  the  input  of  the  plant 
through  a  weighted  input  complementary  sensitivity,  and  recover  tracking  performance 
through  weighted  output  sensitivity.  These  must  be  done  while  keeping  the  robustness  at 
the  output  of  the  plant. 


Figure  64  Time  Responses  due  to  a  step  ac  and  0C 
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Now  a  multiplicative  uncertainty  block  at  the  input  of  the  plant  is  assumed.  The 
multiplicative  uncertainty  represents: 

1.  Uncertainty  in  the  canard  and  the  elevon  actuators 

2.  Uncertainty  in  the  force  and  moments  generated  on  the  aircraft,  due  to  specific 
deflection  of  the  canard  and  elevon 

3.  Uncertainty  in  the  linear  and  angular  accelerations  produced  by  the 
aerodynamically  generated  forces  and  moments 

4.  Others  forms  of  uncertainty  that  are  less  well  understood.  [BDGPS91] 

Figure  6-3  shows  the  block  diagram.  The  dotted  block  represents  the  true  plant.  The 
transfer  function  A(s)  is  assumed  to  be  stable,  unknown,  and  with  an  infinity  norm  less 
than  one  ( II  A(s)  ll„  <  1).  The  objective  for  this  design  is  to  meet  Robust  Stability  at  the 

input  of  the  plant. 


Figure  6-5  Block  Diagram  of  multiplicative  uncertainty  at  the  plant  input 
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Robust  Stability  at  the  input  of  the  plant  is  met  if 

il  A(s)Tedi  (s)  IL  <  1 


by  the  small  gain  theorem.  Thus,  if  we  satisfy 

II  ll„<  1/11  A(s)  IL 


we  have  Robust  Stability  at  the  input  of  the  plant,  or  since  II  A(s)  ll„  <  1, 

HT-iIL<1  (6-3) 

where  Ta,t  is  the  weighted  closed-loop  transfer  function  from  d,  to  e,.  Now  we  have 


Tea,  =  W*,(.)  Tj(s)  (6-4) 


Wdd(s),  called  the  uncertainty  weight,  represents  a  stable  transfer  function  of  the  form 

Wdel(s)  =  wdd(s)*I(2) 


_  50  *  (s  +  100) 
”  (s  +  10000) 


(6-5) 


The  weighting  function  is  used  to  normalize  the  size  of  the  unknown  perturbation  A.  At 
any  frequency  0),  the  value  of  Iw^co)!  can  be  interpreted  as  the  percentage  of  uncertainty 
in  the  model  at  that  frequency.  The  particular  uncertainty  weight  chosen  for  this  problem 
indicates  that  at  low  frequencies,  there  is  potentially  a  50%  modeling  error,  and  at  a 
frequency  of  173  rad/sec,  the  uncertainty  in  the  model  is  up  to  100%  [BDGPS91].  Figure 
6-6  shows  the  Bode  magnitude  plot  of  wdel(s).  The  H„  design  plant  is 


> 

J 

B  ' 

u-l 

= 

c*. 

D* 
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t 

Q 

(6-6) 
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Figure  6-6  Magnitude  of  Multiplicative  Uncertainty  Weighting  Function  (dB) 


and  the  state  space  matrices  are 


[■;]  ■  ft  -I:]  *  ft}*  *  £)-’ 

[ej  =  [  0  C,(  X*  j  +  [  0  (d,]  +  [Ddeiju]  (6-7) 

[y]  =[c8  °(2]  +1°W  +  [°Iu] 

Since  is  not  full  rank,  this  is  a  singular  H,.  design.  The  setup  for  the  mixed 

problem  is:  Find  a  stabilizing  compensator  that  achieves 

jUT~L  s“bi«,R>  HT-lsr,  (6-8) 


K  stabilizing 
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The  performance  index  in  the  numerical  method  is 

(6-9) 

Two  controller  orders  are  presented.  The  first  one  is  a  fourth  order  controller;  it 
represents  the  order  of  the  H2  design.  The  second  one  is  a  sixth  order  controller  and 
represents  the  order  of  the  full  information  Hj/H^  plant.  Therefore,  the  6th  order 
was  found  by  wrapping  the  weights  of  the  H2  design  and  the  weights  of  the  H_  design  into 
one  system.  The  starting  controller  is  (4th  or  6th  order).  Table  6-2  shows  the 
results. 


Table  6-2  Input  Complementary  Sensitivity  Design  lITjJk  and  ||  W^TJL. 


4th  order  Controller 


6th  order  Controller 


II  Tw  ||2 

II  »'<  jdT,  L 

II  T™  |b 

II  wWTi  ||„ 

4.6970 

37.0491 

4.6970 

37.0491 

4.6970 

34.1981 

4.6970 

34.1981 

63683 

2.1831 

6.5190 

2.0029 

6.4369 

1.8485 

6.3395 

1.4366 

63611 

1.4709 

6.3674 

1.2844 

6.7332 

1.0585 

6.6457 

0.8946 

6.9970 

0.6413 

6.8571 

0.7113 

*  7.0968 

03497 

0.5151 

**  7.7440 

03099 

**  7.6167 

03095 

Notice  from  Table  6-2  that  Robust  Stability  is  met  for  both  controllers,  but  there  is  a 
degradation  in  the  Hj  performance.  Figure  6-7  shows  the  a  vs.  7,  curve  for  the  4th  order 
controller,  which  is  virtually  identical  to  that  of  the  6th  order  controller.  The  vector  gain 
and  phase  margins  at  the  input  and  output  of  the  plant  are  shown  in  Table  6-3  for  the  last 
two  controllers  in  both  cases.  As  expected,  the  system  is  very  robust  at  the  input  of  the 
plant;  also  iK  *ice  that  good  margins  at  the  output  were  preserved. 
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Figure  6-7  ||  vs.  ||  T,  |L  (4th  order  controller) 


Table  6-3  Input  Complementary  Sensitivity  Design  VGM^,  and  VPM^, 

controller  (4th  order) 


VGM,  (dB) 

VPMi(dea) 

VGM„  (dB) 

VPM„  (deg) 

*s 

[-5.3269,  16.2768] 

±50.0787 

[-4.5991,  10.4018] 

±40.8564 

■Oi 

[-20.3299,  5.5921] 

±53.7266 

[-11.5388,  4.7866] 

±43.1300 

**s 

[-5.4318,  17.6493] 

±51.5019 

[-4.0590,  7.8659] 

±34.6569 

**T 

[-40.8853,  6.0000] 

±60.0000 

[-8.7230,  4.2634] 

±36.9444 

K^,  controller  (6th  order) 


VGM,  (dB) 

VPMt(deg) 

VGM„(dB) 

*  S 

[-5.3644,  16.7411] 

±50.5846 

[-4.6992,  10.9871] 

±42.0618 

ID 

[-30.0620,  5.8832] 

±57.9334 

[-13.2685,  5.0228] 

±46.0922 

**s 

[-5.3386,  16.4192] 

±50.2367 

[-4.1372,  8.1814] 

±35.5235 

[-34.6236,  5.9396] 

±58.7752 

[-8.3343,  4.1738] 

±35.9332 
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The  H2  design  showed  that  the  system  does  not  perform  tracking  at  all.  Therefore, 
good  performance  will  be  characterized  in  terms  of  the  H„  norm  of  the  output  sensitivity. 

The  output  sensitivity  will  be  a  weighted  sensitivity  function  as  shown  in  Figure  6-8.  The 
exogenous  input  and  output  of  the  plant  are  denoted  as  and  respectively.  The 
weighted  sensitivity  is  our  Nominal  Performance  requirement.  The  weight  for  sensitivity 
represents  an  output  perturbation.  The  transfer  function  is 

Qj/di  =  WpS0  ;  Wp(s)  =  WpCs)*!^ 


wp(s) 


.5*(s  +  3) 
(s  +  0.03) 


(6-10) 


As  in  the  uncertainty  modeling,  the  weighting  function  Wp  is  used  to  normalize 
specifications;  in  this  case,  to  define  performance  by  whether  a  particular  norm  is  less 
than  1.  Nominal  Performance  is  assured  when 


^L<» 


or 


(6-11) 


Figure  6-8  H_  design  for  Nominal  Performance 
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and  since  wp  is  a  SISO  transfer  function,  the  maximum  singular  value  plot  of  the  output 


sensitivity  transfer  function  must  lie  below  the  plot  of 


at  every  frequency.  That  is,  if 


|wp(/  +  GK)'1!  <  1,  then  at  all  frequencies,  ||(/  +  GK)'1(yo))j|^  <  |l  /  wp(y<D)|.  The 


inverse  of  the  weight  wp  is  shown  in  Figure  6-9.  This  sensitivity  weight  indicates  that,  at 
low  frequencies,  the  closed-loop  system  should  reject  output  disturbances  by  a  factor  of 
50-1  [BDGPS91].  The  closed-loop  system  should  perform  better  than  the  open-loop  for 
frequencies  up  to  1.73  (rad/sec),  and  for  higher  frequencies,  the  closed-loop  performance 
should  degrade  gracefully,  always  lying  underneath  the  inverse  of  the  weight  wp  as  shown 
in  Figure  6-9  [BDGPS91]. 


Figure  6-9  Inverse  of  Performance  Weighting  Function,  wp 
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The  design  plant  is  P„2,  given  by 


ft 

1 

< 

Bn 

U-2 

p 

II 

Dv> 

D,u 

CjU 

c,, 

D*. 

Dy«. 

and  the  state  space  matrices  are 


(6-12) 


[e2]  —  [  DpCg 


+  [0M 


[y]  =t  c. 


+  IHOJ  +  [Olu] 


(6-13) 


Sine*  Dj  D  u  is  not  full  rank,  this  is  a  singular  H„  design.  The  mixed  Hj/H^  problem  is 
now  the  solution  to  a  regular  H2  problem  subject  to  a  singular  design.  The  setup  for 
the  mixed  H2/Hoo2  problem  is:  Find  a  stabilizing  compensatin'  that  achieves 


inf  |T 

K  stabilizing*  2 


*w|l2 


subject  to  |TedJ|_<Y2 


(6-14) 


The  performance  index  in  the  numerical  method  is 


Two  controller  orders  are  presented,  4th  and  6th  order.  The  starting  controller  is 
(4th  or  6th  order).  These  are  the  same  as  generated  in  the  previous  section.  Table  6-4 
shows  the  results. 
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Table  6-4  Output  Sensitivity  Design  [|T„||2  and  ||  WpS0|L 


4th  order  < 

Controller 

6th  order  < 

Controller 

HT1J|a 

■D3E39W 

IITJIj 

■EH9 

4.697 

52.1574 

4.697 

52.1574 

4.6994 

44.1814 

4.6999 

402441 

4.7007 

36.186 

4.701 

36.1963 

4.7027 

28.2262 

4.7012 

28.1995 

4.7062 

20.2295 

4.7063 

20.2322 

4.7129 

12.2376 

4.7134 

12.1711 

4.7498 

4.2903 

5.1424 

5.3273 

5.1901 

1.0049 

5.1647 

1.4016 

5.2453 

0.9647 

5.2568 

1.2114 

*  5.2584 

0.9467 

1.1291 

1 _ 

*  *  5.6429 

0.8311 

Notice  from  Table  6-4  that  Nominal  Performance  is  met  for  both  controller  orders.  The 
degradation  in  the  H2  performance  is  not  too  large.  Figure  6-10  shows  the  a  vs.  y2  curve. 
The  vector  gain  and  phase  margins  at  the  input  and  output  of  the  plant  are  shown  in  Table 
6-5  for  the  last  controllers  in  both  cases.  As  expected,  the  system  is  not  robust  at  the 
input  of  the  plant.  Also,  notice  that  the  margins  at  the  output  were  reduced.  This 
represents  the  trade  off  between  performance  and  robustness. 


Figure  6-10  II ||2  vs.  ||  Wp  S0  |L 
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Table  6-5  Output  Sensitivity  Design  VGM^  and  VPM^ 
controller  (4th  order) 


VGMi  (dB) 

VPMi  (deg) 

VGMtt(dB) 

VPM.(des) 

*s 

[-0.9002,  1.0044] 

±  6.2597 

[-3.6380,  6.3789] 

±30.1522 

msa 

[-1.0015,  0.8979] 

±6.2430 

[-5.5480,  3.3584] 

±27.3035 

K,,^  controller  (6  th  order) 


VGM,  (dB) 

VPMi(deg) 

VGM„(dB) 

VPM„  (deg) 

*S 

[-0.9438,  1.0590] 

±6.5802 

[-3.7318,  6.6827] 

±31.1321 

na 

[-1.0560,  0.9415] 

±6.5629 

[-5.7776,  3.4393] 

±28.1164 

**s 

[  3.9341,  1.0468] 

±6.5089 

[-4.1272,  8.1405] 

±35.4130 

**T 

[-1.0437,  0.9316] 

±6.4904 

[-7.4109,  3.9398] 

±33.3542 

Figures  6-11  and  6-12  show  the  time  responses  due  to  a  step  angle  of  attack  command 
and  a  step  pitch  angle  command  (noise  included  in  simulation)  respectively,  for  the  H2 
design  and  the  mixed  designs.  Notice  that  the  design  does  not  provide 

performance  at  all,  which  is  more  evident  in  the  pitch  angle.  The  mixed  controller  shows 
an  improvement  in  the  tracking  performance,  and  the  degradation  in  the  noise  rejection  is 
not  considerable.  The  6th  order  controller  did  not  significantly  improve  the  infinity  norm 
of  the  robustness  and  performance  objectives  nor  the  two  norm  of  T^;  therefore,  a  fourth 
order  controller  seems  to  be  the  best  solution  of  the  two  for  this  mixed  If/H^  design 

problem. 
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This  section  will  examine  the  trade  off  between  weighted  output  sensitivity  (T^  and 
weighted  input  complementary  sensitivity  (Tedl).  Since  the  H2  design  was  shown  to  have 
good  margins  at  the  output  of  the  plant,  let's  examine  the  weighted  output  complementary 
sensitivity.  In  other  words,  the  objective  is  to  observe  how  the  robustness  at  the  output  of 
the  plant  is  affected  by  T^i  and  T^.  The  weight  for  the  output  complementary  sensitivity 
is  chosen  to  be  the  same  as  the  weight  for  input  complementary  sensitivity,  but  in  this  case 
it  is  assumed  to  be  a  fictitious  uncertainty  block  at  the  output  of  the  plant 
[G(I  +  ^dti^fictaious  )!•  Tim  means  that  the  true  plant  is  represented  by  an  input 

multiplicative  perturbation  [(I  +  W^AjG]  and  the  performance  objective  is  weignted 


output  sensitivity.  The  block  diagram  for  output  complementary  sensitivity  is  shown  m 
Figure  6-13. 


u 

Figure  6-13  Weighted  Output  Complementary  Sensitivity 
Block  Diagram 
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Thus, 


T-f=WiliwT#w  (6-16) 

Figure  6-14  shows  how  the  minimization  of  weighted  input  complementary  sensitivity 
(Mixed  )  affects  the  infinity  norm  of  weighted  output  sensitivity  (Mixed  H2/H„2  ) 

and  how  the  minimization  of  weighted  output  sensitivity  affects  the  infinity  norm  of 
weighted  input  complementary  sensitivity.  The  curve  Hj/H^  shows  that  the  minimization 

of  weighted  output  sensitivity  drives  the  infinity  norm  of  weighted  input  complementary 
sensitivity  to  smaller  values.  The  curve  shows  that  the  minimization  of  the 

weighted  input  complementary  sensitivity  starts  to  minimize  the  infinity  norm  of  weighted 
output  sensitivity  also,  but  when  IM_  reaches  small  values,  it  causes  an  increase  in 

|WpS0|  .  Figure  6-15  shows  how  the  minimization  of  weighted  input  complementary 

sensitivity  does  not  affect  the  weighted  output  complementary  sensitivity.  This  means  that 
this  design  does  not  affect  the  robustness  at  the  output.  The  values  of  the  infinity  norms 
for  the  different  transfer  functions  are  in  Appendix  A,  Sectio  %  A.  1  and  A.2. 


Figure  6-14  vs.  J\VpS0|L  for  Hj/H^  design  and  H2/H-2  design  (K  4th) 
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Figure  6-15  |W*,T#|L  vs.  JW^TjL  for  design  (K  4th) 

Figure  6-16  shows  how  the  minimization  of  weighted  output  sensitivity  affects  the 
weighted  output  complementary  sensitivity.  Notice  how  the  robustness  at  the  output  of 
the  plant  starts  to  decrease  as  the  system  gets  more  performance.  The  conclusion  is  that 
when  the  weighted  output  sensitivity  is  reduced,  it  drives  the  weighted  complementary 
sensitivity  to  higher  values. 


Figure  6-16  vs.  Ks-L  for  Hj/H^  design  (K  4th) 
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The  setup  for  this  mixed  Hj/H^  problem  is:  Find  an  stabilizing  compensator  that 


achieves 


*.2L,IT»IL  subjectto 


fllW^T, 


wPs„ 


L  -Yi 
~  Y2 


(6-17) 


where  both  constraints  will  be  treated  as  equality  constraints.  The  performance  index  for 
die  numerical  method  is 

■/,  =K&  +^i(l'VddTiJ_  -  Y,)2  +  ^KS4. -Y,)2  (6-18) 


The  state  space  matrices  are  equation  (6-2)  for  the  Hj  part,  equation  (6-7)  for  the 
weighted  input  complementary  sensitivity,  and  equation  (6-13)  for  the  weighted  output 
sensitivity.  This  design  will  map  the  boundary  between  |WdelTi|L  |WPSo|L  When 

both  infinity  norms  are  close  to  the  optimal  values,  respectively.  This  will  be  done  by 
minimizing  both  constraints  using  the  direct  method  in  different  directions  as  shown  in 
Figure  6-17.  Two  cases  are  defined.  Case  1  tries  to  reduce  as  much  as  possible  the 
infinity  norm  of  the  weighted  output  sensitivity  while  holding  the  infinity  norm  of  the 
weighted  input  complementary  sensitivity  less  than  one.  In  other  words,  the  first  case  tries 
to  get  the  best  level  of  performance  that  meets  the  robustness  requirement.  Case  2  tries  to 
reduce  both  infinity  norms  as  much  as  possible,  which  means  that  it  is  desired  to  get  the 
best  performance  and  the  best  robustness.  The  infinity  norm  of  the  weighted  output 
complementary  sensitivity  will  also  be  calculated  in  order  to  observe  the  trade  off  between 
this  design  and  the  robustness  at  the  output  of  the  plant 
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Figure  6-17  Objectives  of  the  mixed  problem  with  two  H„  constraints 
Again  the  starting  controller  is  the  optimal  Hj  controller,  which  is  the  order  of  the  Hj 
part  only  (fourth  order);  then  the  controller  order  is  increased  to  6th  (explained  before) 
and  8th  order  (computed  by  wrapping  P2,  P„,,  and  P„2  into  one  system,  P).  Therefore, 

fourth,  sixth,  and  eight  order  mixed  controllers  will  be  generated.  The  method  used  in  the 
numerical  technique  is  the  direct  method.  Table  6-6  shows  part  of  the  results  (see 
Appendix  A,  Section  A.3  for  more). 


lITJk 

IIW^T.IL 

4.6970 

37.0491 

52.1574 


Table  6-6  Mixed  H/H.  with  two  H.  Constraints:  | JTJIj ,  ||  W^TjL,  and  ||  W_S0  |L 


IIW^TJL 


0.6853 


Fourth  order  controller 


0.9938  0.9234 


0.8172  0.9110 


Sixth  order  controller 


0.9978  0.9574 


Case  1 


Case  2 


5.6883 


6.1322 


1.0159 


0.9971 


Case  1 


Case  2 


5.7311 


6.0274 


0.8519 


0.9579 


0.9826 


0.9828 


Eight  order  controller 


Case  1 

5.9228 

0.9220 

0.6845 

0.7560 

** 

5.8145 

0.7438 

0.7644 

0.7030 

Case  2 

5.6926 

0.6109 

0.7270 

0.6707 
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Figure  6-18  Direct  Method  versus  Single  H„  constraint  designs  (HIMAT ) 
Notice  that  the  6th  order  controller  achieves  similar  results  to  the  4th  order,  and  the 
robustness  at  the  output  of  the  plant  is  maintained  for  both.  The  big  Difference  is  the  8th 
cider  controller,  which  achieves  smaller  infinity  norms  than  the  4th  and  6th  order 
controller,  and  also  keeps  the  level  of  robustness  at  the  output  (the  **  8th  order  controller 
will  be  explained  in  Section  6.7).  Figure  6-18  shows  how  the  direct  method  with  multiple 
constraints  goes  directly  to  the  minimum  values  of  the  constraints.  This  figure  is 

similar  for  the  4th,  6th  and  8th  order  controllers  using  die  direct  method.  Notice  that  from 
a  control  point  of  view,  we  are  not  interested  in  values  of  the  HL  constraints  above  one, 

since  they  do  not  meet  out  requirements.  The  following  figures  will  show  the  area  of 
interest  (below  one)  for  the  H_  constraints.  Figure  6-19  shows  the  trade  off  between 

W^Tj  and  WpSo  using  a  4th  order  controller.  The  constraint  boundary  between  weighted 
input  complementary  sensitivity  and  weighted  output  sensitivity  is  defined.  This  boundary 
shows  that  Case  2  is  better  than  Case  1,  since  it  obtains  better  values  of  performance  and 
robustness.  The  6th  order  controller  shows  similar  results  to  the  4th  order  controller. 
Figure  6-20  shows  the  trade  off  between  W^Tj  and  WpSo  using  an  8th  order  controller. 
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Notice  that  this  order  of  controller  improves  the  values  of  the  infinity  norms,  and  the  non- 
feasible  region  is  also  "smaller".  Appendix  A,  Section  A.3  shows  the  complete  table  of 
results. 


Figure  6-20  Trade  off  between  weighted  T,  and  S0  (8th  order  controller) 


Figure  6-21  WJs)  T,(s)  Maximum  Singular  values  (K^  8th  order) 


The  maximum  singular  values  of  WJsJTjCs)  are  shown  in  Figure  6-21.  The  y  axis 
represents  some  of  the  different  8th  order  controllers;  it  starts  at  which  corresponds 
to  y  =  1,  and  moves  into  the  mixed  controllers  (from  y  =  2  through  y  =  13).  Notice  how 
the  maximum  singular  values  are  minimized  until  Robust  Stability  is  met.  Figure  6-22 
shows  the  maximum  singular  values  of  Wp(s)S0(s).  The  same  axis  orientation  is  kept. 
This  new  technique  drives  both  infinity  constraints  to  the  required  level  of  Nominal 
Performance  and  Robust  Stability,  with  the  two  norm  also  being  minimized.  It  is 
important  to  mention  again  that  both  designs  are  singular  problems.  Table  6-7  shows 

the  VGM  and  VPM  for  each  controller  using  the  final  controller.  Notice  that  each 
controller  tries  to  recover  the  VGM  and  VPM  at  the  input  of  the  plant.  Now  the  trade  off 
between  input  and  output  margins  is  more  evident,  since  a  weighted  output  sensitivity 
output  is  the  driver  for  margins  at  the  output  of  the  plant. 
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1— >13  K2opl~»  Kmix 


Frequency  (redAec)  (1 0*xpi-{-3,3]) 

Figure  6-22  Wp(s)S0(s)  Maximum  Singular  values  (K^  8th  order) 


Table  6-7  Mixed  Hj/Hw  with  two  H_  Constraints  VGM  and  VPM 


_ £ 

*  controller  ( 

4th  order) 

VGM,  (dB) 

VGM„(dB) 

SI 

[-3.9103,  7.3025] 

±33.0344 

[-3.7631,  6.7876] 

±31.4629 

T1 

[-6.3912,  3.6419] 

±30.1926 

[-6.4753,  3.6682] 

±30.4666 

S2 

[-4.3418,  9.0814] 

±37.8401 

[-3.8033,  6.9242] 

±31.8883 

T2 

[-8.7051,  4.2594] 

±36.8988 

[-6.4633,  3.6645] 

±30.4278 

_ £ 

ri,  controller  (6th  order) 

SI 

±29.9668 

[-3.6540,  6.4297] 

±30.3183 

T1 

[-6.2052,  3.5825] 

±29.5776 

[-7.2142,  3.8858] 

±32.7709 

S2 

[-4.1892,  8.3994] 

±36.1056 

[-3.6533,  6.4275] 

±30.3112 

T2 

[-8.1393,  4.1269] 

±35.4095 

[-7.1817,  3.8768] 

±32.6735 

_ K, 

controller  (1 

8th  order) 

SI 

[-3.6560,  6.4359] 

±  30.3386 

[-4.7985,  11.6176] 

±43.2774 

T1 

[-6.8178,  3.7721] 

±31.5575 

[-9.8992,  4.5066] 

±  39.7588 

s** 

[-4.5405,  10.0795] 

±40.1593 

[-4.4862,  9.7928] 

±39.5186 

p»# 

[-8.7010,  4.2584] 

±36.8883 

[-7.3548,  3.9246] 

±33.1891 

S2 

[-4.7691,  11.4256] 

±42.9160 

[-4.7003,  i 0.9939] 

±42.0753 

T2 

[-14.9428,  5.2062] 

±48.4723 

[-12.1456,  4.8756] 

±44.2334 

(1=  case  1, 2=  case  2) 
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This  table  suggests  the  use  of  a  higher  order  controller  in  order  to  avoid  a  degradation  of 
margins  at  the  output  of  the  plant.  The  next  section  will  compare  the  (**)  eight  order 
mixed  controller  with  the  controller  obtained  using  p  -synthesis,  as  well  as  an  optimal 

controller. 

Miredlfc/Hoo  with  Multiple  Singular  H~>  Constraints 

and  Robust  Performance 
fi-Svnthesis  versus  Mixed  H?/H^-  Synthesis 

This  section  presents  a  comparison  between  the  p  controller  obtained  from  running 
the  HIMAT  demo  in  MATLAB,  the  FL  optimal  controller  for  the  mixed  sensitivity 
problem,  and  the  mixed  H2/H„  controller.  In  this  section,  the  value  of  the  upper  bound  on 

p.(a(DMD_1 ))  will  simply  be  called  p  for  convenience.  A  fictitious  block  Ap  is  created 
to  include  the  performance  requirement,  as  shown  in  Figure  6-23. 


Figure  6-23  Block  diagram  for  Robust  Performance 
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Although  the  mixed  problem  does  not  directly  address  Robust  Performance,  the  trade  off 
among  the  infinity  norms  of  the  diagonal  and  cross  terms  of  M  will  shape  the  maximum 
singular  values  of  M,  and  therefore  could  affect  p.  This  is  shown  when  a  p  analysis  is 
done  for  twenty-six  8th  order  mixed  controllers  from  Section  6.6  (these  26  controllers  are 
only  a  subset  of  all  8th  mixed  controllers;  see  Appendix  A,  Section  A.3),  as  shown  in 
Figure  6-24.  is  defined  as  the  upper  bound  on  p(M).  The  x  axis  corresponds  the 

different  controllers,  starting  with  and  moving  along  the  mixed  controllers  as  the 
infinity  norms  are  minimized.  Notice  that  the  upper  bound  of  p  is  minimized.  Now, 
looking  at  the  three  8th  order  mixed  controllers  from  Table  6-6,  we  see  that  the  case  1  and 
case  2  controllers  have  a  IIMII^  value  around  1.7,  while  the  controller  (**)  has  IIMII^  = 

1.3404,  as  shown  in  Figure  6-25.  This  shows  that  different  combinations  of  infinity  norms 
of  the  cross  and  diagonal  terms  will  result  in  a  differing  values  of  IIMII^.  This  is  the  reason 

for  using  the  D  scaling  (see  Chapter  2),  since  it  shapes  the  maximum  singular  values  of  M 
in  order  to  reduce  IIMII^.  The  mixed  problem  improves  IIMIl^,  but  it  does  not  directly 

address  Robust  Performance,  since  it  does  not  exploit  the  frequency  information. 
Consider  the  controller  (**)  as  the  best  controller  in  terms  of  Robust  Performance.  Figure 
6-26  shows  the  p  bounds  for  the  p-synthesis  design,  H..  optimal  control,  and  the  (**) 
mixed  controller.  The  mixed  controller  reduces  IIMIIH  more  than  the  H„  optimal  controller 
does,  because  it  directly  addresses  Robust  Stability  and  Nominal  Performance. 
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Figure  6-24  Mixed  H^/H.  controllers  and  the  effect  on  ||  M  1^ 


Figure  6-25  m.(M)  of  the  8th  order  Mixed  controllers  from  Table  6-6 


Figure  6-26  fi(M)  of  the  ^-synthesis,  H„  ^  and  mixed  H2/H_ControUer 

control  design 


Although  the  mixed  Ha/H^controller  does  not  pass  the  test  for  Robust  Performance,  the 
order  of  this  controller  is  much  smaller  than  the  K(s)  obtained  from  p-syn thesis.  Also, 
noise  rejection  is  better  with  the  mixed  controller  than  the  other  two.  Table  6-8 
summarizes  the  results,  including  a  4th  order  mixed  controller  as  well. 


Table  6-8  |i  Analysis  for  the  HIMAT  example 


Robust 

Stability 

Nominal 

Performance 

Robust 

Performance 

invji* 

order  of 
controller 

pass 

pass 

0.9803 

1851034 

20 

Mixed  Hj/H^ 

pass 

pass 

1.3404 

5.81450 

8 

*H„  optimal 
controller 

pass 

pies 

1.6230 

oo 

7 

Mixed  ItyH.. 

pass 

5.4040 

6.1320 

4 

*  the  optimal  H^,  controller  is  for  mixed  sensitivity  and  complementary  sensitivity 
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Figures  6-27  and  6-28  show  the  magnitude  of  the  weighted  sensitivity  function  and 
weighted  complementary  sensitivity  function  for  the  different  control  designs,  respectively. 
Notice  how  the  infinity  norm  for  both  functions  are  improved  for  the  mixed  problem 
compared  to  the  problem.  Figures  6-29  and  6-30  represent  the  magnitude  plot  versus 

frequency  for  the  actual  (unweighted)  sensitivity  and  complementary  sensitivity  functions. 
Figure  6-29  shows  that  the  mixed  Hj/H^  controller  tries  to  decrease  the  magnitude  of  the 

sensitivity  function  at  low  frequencies,  compared  to  the  p  controller.  This  is  probably  due 
to  the  fact  that  the  mixed  controller  is  trying  to  reject  the  low  frequency  wind  disturbance 
which  p  does  not  account  for.  In  Figure  6-30,  the  mixed  Hj/H^  rolls  off  faster  than  p 
does,  and  peaks  earlier  than  H_. 


Figure  6-27  Magnitude  of  Wp(s)S0(s)  for  p-Synthesis,  H„,  and  the  Mixed  Hj/H. 
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Figure  6-28  Magnitude  of  WM(s)T,(s)  for  ^-Synthesis,  H..,  and  Mixed  Hj/H.. 


Figure  6-29  Magnitude  of  Sensitivity  Function  for  ^-Synthesis,  Hm, 

and  Mixed  H2/H„ 
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Figure  6-30  Magnitude  of  Complementary  Sensitivity  for  p-Synthesis, 

H„,  and  Mixed 


The  time  responses  for  a  step  angle  of  attack  and  pitch  angle  are  shown  in  Figures  6-31 
and  6-32.  Notice  how  the  wind  disturbance  and  noises  affect  the  p-syn  thesis  controller 
and  the  controller.  The  mixed  Ha/KL  controller  has  good  noise  rejection  and  good 

tracking  as  well  The  mixed  controller  fails  the  robust  performance  test  with  its  upper 
bound  of  1.3404;  however,  this  means  that  in  order  to  pass  this  test,  the  robustness  and 
nominal  performance  requirements  must  be  relaxed  by  only  a  factor  of  1/1 .3404. 
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<1.8  Mixed  design  witli  three  li»>  constraints 

The  objective  of  this  design  is  to  minimize  three  constraints.  The  first  and  second 

constraints  are  the  previous  Robust  Stability  and  Nominal  Performance  requirements. 
These  H_  constraints  were  examined  individually  and  simultaneously  in  the  previous 

sections.  Now  a  weighted  output  complementary  sensitivity  constraint  is  introduced,  as  a 
third  constraint.  The  desire  is  to  drive  its  infinity  norm  to  less  than  one.  The 
exogenous  input  is  d3  and  the  exogenous  output  is  for  this  new  H„  constraint,  as  shown 
in  Figure  6-33.  The  dotted  block  represents  our  true  plant  The  transfer  function  Af(s)  is 

a  fictitious  unstructured  block  and  is  assumed  to  be  stable,  unknown,  and  such  that  its 
infinity  norm  is  less  than  one  ( II  A£s)  ll„  <  1).  The  objective  for  this  design  is  to  meet 

II IL<1  ;  Robust  Stability  at  input 
H  Tod2IU<  1  ;  Nominal  Performance  at  output  (6-19) 

«T*L<  1;  Robust  Stability  at  output 
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The  state  space  for  the  H^3  design  plant  is 


+ GN + [o]u] 

[ej]  =  [  D«C.  &-{*]  +  [Old,]  +  [  0  Iu]  (6-21) 

[y]  =[  c.  o{*]  +  HM  *  [«M 

* 

is  not  lull  rank;  therefore,  this  is  a  singular  H„  problem.  The  setup  for  this 
mixed  problem  is:  Find  a  stabilizing  compensator  that  achieves 


(jW^Ti  M  <  y, 

subject  to  ■  jw,  S.|_  <iy,  (6-22) 

llWMT„LsY, 

Again,  the  three  H_  constraints  will  be  treated  as  equality  constraints.  Thus,  the 
performance  index  for  the  numerical  method  is 

A  =|r_£ +A(K,t,|_ -y,)2 +x2(|wps0|_ -y2)2  +  -  y,)2 

(6-23) 

The  state  space  matrices  are  equation  (6-2)  for  the  Hj  part,  equation  (6-7)  for  weighted 
input  complementary  sensitivity,  equation  (6-13)  for  the  weighted  output  sensitivity 
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output  and  equation  (6-21)  for  weighted  output  complementary  sensitivity  for  the  H„ 
parts. 

Again,  the  starting  controller  is  the  optimal  Hj  controller.  The  design  is  performed  to 
obtain  fourth,  sixth,  and  eight  order  controllers.  The  method  used  in  the  numerical 
technique  is  the  direct  method.  Table  6-9  shows  part  of  the  results  (see  Appendix  A, 
Section  A.4  for  foil  results).  Comparing  Table  6-9  with  Table  6-6  (two  H„,  constraints), 

we  can  see  that  there  is  no  improvement  for  the  4th  and  6th  order  controllers.  It  seems 
that  the  third  constraint  is  not  dominant,  and  the  trade  off  among  the  H„  constraints  is 

only  related  to  the  weighted  input  complementary  sensitivity  and  weighted  output 
sensitivity.  The  8th  order  controller  does  show  some  improvement  in  performance,  while 
keeping  a  good  level  of  robustness  at  the  input  and  output  of  the  plant  Notice  that  if  an 
controller  is  found  for  an  unmixed  problem,  the  order  would  be  10.  A  tenth  order 

mixed  controller  was  not  computed  due  to  time  constraints. 


Table  6-9  Mixed  Hj/H„  with  three  H„  Constraints:  IfT^Hj,  ||  WMTJL, 


||W;S,|Uand||WwT0|L 


1  IT.I, 

IIW*,T,|L 

II  WwT#  |L 

IC,„.  Controller 

4.6970 

tHEZdHI 

52.1574 

0.6853 

j 

Fourth  order  controller 

I  6.5415 

!  0.7584 

0.9283 

0.9975 

j 

Sixth  order  controller 

!  6.5775 

0.8029 

0.9416 

0.9894 

j 

Eight  order  controller 

5.2653 

i  0.9189 

0.6865 

0.5956 

j 

Table  6-10  shows  the  VGM  and  VPM  at  the  input  and  output  of  the  plant.  The  time 
responses  for  a  step  angle  of  attack  (a)  and  pitch  angle  (0)  command  are  shown  in  Figures 
6-34  and  6-33  for  the  controllers  from  Table  6-9.  Notice  that  the  4th  and  6th  order 
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controllers  have  a  considerable  overshoot  Although  the  settling  time  is  almost  the  same 
for  all  three,  the  8th  order  controller  has  a  better  response. 


Table  6-10  Mixed  ItyH..  with  three  H_  Constraints:  VGM  and  VPM 

controller  (4th  order 


1  VGM.  (dB) 

VPM,  (dea)  VGM„  (dB) 

VPM.  (deg) 

S 

[-4.2084,  8.4821] 

±36.3227 

[-3.7627,  6.7862] 

±31.4586 

T 

[-9.9643,  4.5189] 

±39.9042 

[-6.3492,  3.6287] 

±30.0548 

[-4.2646,  8.7285] 


[-8.8461,  4.2907] 


[-3.7919,  6.8851] 


[-6.8231,  3.7736] 


controller  (6th  order 


±  36.9585 


±  37.2559 


[-3.7126,  6.6194] 


[-6.4819,  3.6703] 


controller  (8th  order 


±31.7673 


±31.5740 


[-4.9276,  12.5242] 


[-16.1158,  5.3134] 


±30.9307 


±30.4881 


±44.8856 


±49.8972 


M  Summary 

If  the  selection  of  the  best  controller  has  to  be  made,  the  question  by  itself  is  too 
complex,  since  many  factors  have  to  be  considered.  Table  6-11  summarizes  the  most 
important  factors  to  consider  in  the  selection  of  the  controller.  Consider  the  following 
factors  (x  in  Table  6-1 1  if  passed): 

1.  Satisfy  Robust  Stability  (RS)  at  the  input  and  output  of  the  plant,  and  satisfy 
Nominal  Performance  (NP). 

2.  Satisfy  RS  and  NP  with  high  noise  rejection 

3.  Satisfy  Robust  Performance  <  1 

4.  Relax  the  Robust  Performance  <1/1 .3404 

5.  Relax  the  Robust  Performance  <  1/5.4020 

6.  Low  overshoot 
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7.  Low  order  controller  that  satisfies  RS  and  NP 


Table  6-11  (HIMAT)  Selection  of  Controllers 


Factor 

Controller 

1 

2 

3 

1 

5 

6 

1 

Mixed  4th  (two  H..  Const) 

H 

■ 

■ 

■ 

■ 

■ 

H 

Mixed  8th  (two  H_  Const.) 

H 

■ 

H 

H 

■ 

■ 

Mixed  4th  (three  Const.) 

H 

■ 

■ 

■ 

H 

■ 

H 

Mixed  8th  (three  H_  Const) 

H 

H 

■ 

■ 

H 

H 

■ 

H-(8th) 

H 

■ 

■ 

■ 

H 

■ 

■ 

|X  (20th) 

H 

■ 

H 

H 

H 

■ 

■ 

These  are  only  a  few  factors  to  be  considered.  The  selection  will  depend  on  how  well  we 
know  the  plant,  and  what  factors  are  more  important  than  others.  The  order  of  the 
controller  is  very  important  when  it  has  to  be  implemented;  therefore,  the  mixed  Hj/H^ 
control  problem  with  multiple  constraints  will  have  great  importance  in  low  order 
controller  design. 
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VII.  Summary.  Conclusions  and  Recommendations 

Summary 

The  mam  objective  of  this  thesis  was  to  investigate  the  mixed  Hj/H..  with  multiple  H„ 
constraints  control  problem.  A  SISO  and  a  MIMO  problem  were  solved.  Successful 
strategies  for  obtaining  measures  of  performance  and  robustness  were  presented  in  most 
designs. 

The  initial  chapter  gave  a  limited  synthesis  history  that  represented  the  motivation  for 
the  use  of  multiple  H„  constraints  in  the  mixed  problem.  Chapter  II  discussed  the  three 
base  methodologies  (lij,  H_,  and  ^-synthesis),  and  a  review  of  the  related  design 
examples.  Chapter  in  discussed  the  mixed  H2/HL  problem  with  a  non-singular  H„ 
constraint  and  with  a  single  singular  constraint  Also,  the  new  numerical  method  was 
explained.  Chapter  TV  developed  the  mixed  H^..  problem  with  multiple  H„  constraints, 
and  discussed  how  the  new  numerical  method  can  be  used  to  solve  this  problem. 

Chapter  V  presented  the  F-16  short  period  approximation  (plus  first  order  servo  and 
Padd  approximation)  SISO  example.  The  SISO  example  represented  an  introduction  to 
this  new  design  technique.  First  an  Hj  design  was  accomplished,  and  then  two  mixed 
problems  with  one  H.  constraint  were  solved.  The  H„.  constraints  were  weighted  input 
sensitivity  and  weighted  input  complementary  sensitivity.  The  trade  off  between  both  H. 
constraints  was  observed.  Finally,  the  mixed  ItyH^  with  multiple  constraints  problem 

was  solved.  Two  methods  were  applied:  the  grid  method  and  the  direct  method.  A 
surface  was  created  using  the  grid  method,  and  it  showed  possible  boundaries  between  H. 
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constraints  when  they  are  close  to  optimal  values.  The  direct  method  was  shown  to  be  a 
better  method,  since  it  permitted  selection  of  the  direction  of  minimization. 

Chapter  VI  presented  the  HIMAT  problem  (MIMO).  Using  the  HIMAT  problem,  the 
following  were  solved: 

•  An  H2  design. 

•  Two  mixed  Hj/H,.  designs,  each  with  one  H„  constraint  The  H_  constraints 

addressed  Robust  Stability  and  Nominal  Performance  independently. 

•  A  mixed  Hj/H..  design  with  two  H„  constraints.  The  constraints  addressed 

Robust  Stability  and  Nominal  Performance,  but  now  the  trade  off  between  them 
was  manipulated. 

•  An  H„,  -synthesis  and  p-synthesis  (using  D-K  iteration),  for  the  augmented  system 
(the  two  H„.  constraints  wrapped  in  one  block)  were  solved. 

•  Finally,  a  three  H„  constraint  problem  was  solved.  The  constraints  addressed 
Robust  Stability  at  the  input  and  output  of  the  plant,  and  Nominal  Performance. 

Different  order  mixed  controllers  were  produced,  with  good  results  in  most  of  them.  A 
p  analysis  was  done  on  the  mixed  controller,  and  this  was  compared  with  the  optimal 

controller  and  the  p-synthesis  controller.  Table  6-11  summarizes  some  of  the  major 
results. 

Conclusions 

The  conclusions  of  this  thesis  are: 

i) .  This  new  numerical  technique  permits  minimization  of  the  two  norm  of  Tw  subject 
to  single  or  multiple  constraints. 

ii) .  The  constraint  can  be  regular  or  singular. 
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iii) .  As  the  order  of  the  controller  was  increased,  better  results  were  obtained. 

iv) .  The  mixed  FtyH..  optimization  with  multiple  singular  H„  constraints  permits  the 

designer  control  over  the  level  of  Robust  Stability  and  Nominal  Performance,  and  also  to 
fix  the  level  of  other  H„  constraints.  This  means  that  the  trade  off  between  design 
requirements  in  terms  of  constraints  can  be  freely  chosen  by  the  designer. 

v) .  The  main  idea  to  include  other  H„  constraints,  that  are  not  specified  as  a  Nominal 
Performance  or  Robust  Stability  requirement,  is  that  the  system  could  meet  Robust 
Stability,  Nominal  Performance  and  even  Robust  Performance  for  a  certain  number  of 
uncertainty  blocks  and  performance  requirements,  yet  fail  when  the  number  of  uncertainty 
blocks  or  performance  requirements  are  increased.  Therefore,  this  technique  can  keep  the 
level  of  Robust  Stability  and  Nominal  Performance  for  the  original  blocks,  and  can  control 
die  level  in  terms  of  H„.  magnitude  for  those  that  the  system  does  not  meet. 


Table  7-1  summarizes  the  improvement  of  this  new  nonconservative  method, 
compared  with  other  control  design  methods. 

Table  7-1  Comparison  of  different  Control  Law  Designs 


H2 

H„ 

P(D-K) 

Mixed 

h2/pl 

Handle  white  Gaussian  noise  (WGN) 

X 

X 

Robust  Stability,  Nominal  Performance 

X 

X 

X 

Robust  Performance 

X 

Trade  off  between  RS  and  NP  freely 

X 

WGN  and  RS,  NP 

X 

Reduced  order  controller 

X 

'-3 


i)  Improve  the  numerical  method,  especially  around  the  knee  of  the  a  vs.  y  curve. 

ii)  Investigate  the  "relationship"  between  the  diagonal  and  cross  terms  in  mixed  H2/H„ 
with  multiple  H„  constraints,  and  how  this  relationship  could  affect  Robust  Performance. 

iii)  A  faster  computer  is  needed  in  order  to  obtain  a  large  number  of  controllers,  so  that 
various  trade-offs  can  be  examined  quickly  and  efficiently. 

iv)  Since  the  numerical  method  runs  with  any  order  controller,  investigate  results  using 
any  order  stabilizing  controller,  including  order  less  than  the  underlying  H2  plant  (if  it 
exists). 

v)  Remove  the  restriction  that  the  H„  constraints  must  be  satisfied  with  equality 
through  the  use  of  constrained  optimization.  This  has  already  be-n  done  [Wal94],  and  the 
results  in  this  thesis  are  being  reworked  using  sequential  quadratic  programming  rather 
than  DFP. 


A.1  Weighted  Input  Complementary  Sensitivity 


4th  order  controller 


«T-r"a 

HW^TjIL 

IIWpS0IL 

4.6970 

37.0427 

52.1864 

0.6856 

4.6970 

34.1942 

51.7626 

0.6602 

4.7031 

28.1678 

51.9762 

0.5881 

4.8223 

22.2256 

44.8161 

0.6134 

4.8698 

19.5766 

41.2103 

0.5564 

5.0489 

16.3526 

33.9038 

0.5556 

5.2281 

13.4063 

23.7194 

0.6010 

5.4864 

7.4497 

15.1682 

0.6236 

6.3685 

2.1831 

21.3140 

0.6466 

6.4369 

1.8485 

22.9333 

0.6571 

6.5611 

1.4709 

20.3251 

0.6554 

6.7332 

1.0585 

26.1706 

0.6613 

6.9970 

0.6413 

34.6669 

0.6784 

7.7440 

0.4956 

48.7810 

0.7959 

7.1843 

0.4918 

51.2928 

0.6487 

7.1759 

0.4920 

51.7376 

0.6478 

rotter  controller 

hw^il 

HWpS0IL 

IIW^TJL 

4.6970 

37.0427 

52.1864 

0.6856 

4.6970 

34.1941 

51.7622 

0.6603 

4.7064 

28.2220 

50.7688 

0.5726 

4.7365 

25.3941 

48.4703 

0.6247 

5.0925 

16.4510 

37.8934 

0.5719 

5.3813 

13.1794 

32.9813 

0.5940 

5.5764 

10.7019 

24.2157 

0.6249 

6.3416 

4.6897 

23.0864 

0.6366 

6.5190 

2.0029 

24.0980 

0.6154 

6.4690 

1.8122 

24.9896 

0.6316 

6.4751 

1.5727 

25.0295 

0.6737 

6.3395 

1.4366 

21.4272 

0.7059 

6.3674 

1.2844 

25.4598 

0.7302 

6.4721 

1.1122 

28.0406 

0.7332 

A-l 


6.6457 

0.8946 

6.8571 

0.7113 

7.0705 

0.5177 

7.2149 

0.5548 

7.6167 

0.5095 

7.0705 

0.5177 

31.8055 

0.7371 

36.4956 

0.7552 

45.8587 

0.6423 

48.3201 

0.7947 

56.5616 

0.8143 

45.8587 

0.6423 

4th  order  controller 


nr- it 

11  WddTi  II- 

IIWpSJL 

IIW^TJL 

4.6970 

37.0427 

52.1864 

0.6856 

4.6990 

25.9277 

48.2140 

0.5763 

4.6994 

24.5831 

44.2054 

0.5439 

4.7002 

23.8210 

40.2266 

0.5406 

4.7007 

22.8746 

36.2079 

0.5410 

4.7019 

22.0211 

32.2226 

0.5440 

4.7027 

21.2916 

28.2454 

0.5509 

4.7042 

20.3073 

24.2368 

0.5628 

4.7062 

18.7839 

20.2418 

0.5714 

4.7129 

16.0872 

12.2381 

0.6305 

4.7327 

15.2283 

8.3599 

0.7095 

4.7498 

10.7949 

4.2906 

0.9023 

5.1901 

5.3058 

1.0481 

1.1369 

5.2453 

4.9608 

1.0108 

1.1111 

5.2584 

4.8214 

0.9905 

1.0972 

5.2617 

4.8043 

0.9887 

1.0960 

6th  onto  grotroUer 

NT^ 

HW^TjlL 

HWpSJL 

IIW^TJL 

4.6970 

37.0427 

52.1864 

0.6856 

4.6984 

25.6282 

48.2167 

0.5744 

4.6999 

23.7559 

40.2441 

0.5407 

4.7010 

22.9563 

36.2181 

0.5412 

4.7012 

21.2190 

28.2184 

0.5525 

4.7063 

18.3781 

20.2468 

0.5713 

4.7097 

17.1423 

16.2076 

0.5955 

4.7134 

15.4557 

12.1765 

0.6364 

5.1647 

11.3612 

1.4016 

1.4939 

5.2568 

7.0627 

1.2114 

1.0669 

5.3035 

4.6762 

1.1291 

1.0641 

5.6429 

4.5730 

0.8311 

0.9451 

A-2 


IIT.II 2 

HW^TVL 

HWpS0IL 

IIW^TJL 

4.6970 

37.0491 

52.1574 

0.6853 

4.9689 

25.9037 

37.9268 

0.7416 

4.9747 

17.5703 

25.2809 

0.5761 

5.0639 

9.0661 

13.1080 

0.6197 

5.1136 

7.2982 

10.1080 

0.6549 

5.2435 

5.6649 

7.3953 

0.6874 

5.3986 

1.5091 

1.0073 

1.0344 

5.5195 

1.2231 

0.9592 

1.0137 

5.5208 

1.2231 

0.9365 

1.0137 

5.5636 

1.1034 

0.9542 

1.0154 

5.5827 

1.0989 

0.9178 

1.0175 

5.5834 

1.0982 

0.9174 

1.0167 

5.5884 

1.0981 

0.9159 

1.0139 

5.6837 

1.0009 

0.9229 

1.0225 

5.6883 

0.9938 

0.9234 

1.0159  . i 

6.0137 

0.8467 

0.9146 

1.0036 

6.0963 

0.8278 

0.9099 

0.9985 

6.1099 

0.8234 

0.9122 

1.0000 

6.1161 

0.8215 

0.9106 

0.9983 

6.1101 

0.8244 

0.9102 

0.9987 

6.1091 

0.8244 

0.9100 

0.9986 

6.1322 

0.8172 

0.9110 

0.9971  . i 

6.2257 

0.8039 

0.9154 

0.9990 

6.8656 

0.6933 

1.0612 

0.9698 

6.9195 

0.6916 

1.0474 

0.9710 

6.9908 

0.6895 

1.0600 

0.9732 

honterw 

«Twii2 

nttplta 

IIW^TjlL 

HWpS0IL 

4.6970 

37.0491 

52.1574 

0.6853 

4.7468 

34.0391 

46.8218 

0.7340 

4.7762 

29.5118 

41.4565 

0.6448 

4.8104 

26.3026 

37.0553 

0.6860 

4.9899 

8.1814 

10.7511 

0.6518 

5.3292 

1.5742 

1.4497 

0.9343 

5.7274 

1.0706 

0.9707 

1.0010 

5.7311 

0.9978 

0.9574 

0.9826  . 

5.9607 

0.8712 

0.9581 

0.9832 

6.0292 

0.8512 

0.9596 

0.9819 

(Case  1) 


(Case  2) 


(Case  1) 
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6.0274 

0.8519 

0.9579 

0.9828 

6.0341 

0.8161 

0.9764 

1.0039 

6.3459 

0.7776 

1.0317 

0.9624 

6.9503 

0.7459 

1.2159 

0.9175 

5.7020 

1.0383 

0.9481 

0.9748 

5.7287 

1.1739 

0.9177 

0.9909 

5.4175 

1.6316 

0.8939 

0.9889 

5.3692 

1.7694 

0.8891 

0.9475 

5.3261 

3.2591 

0.8712 

0.9776 

8th  order  controller 
irr^iij  iiw^iL 

iiwps0il 

IIW^TJI, 

4.6970 

37.0491 

52.1574 

0.6853 

4.7023 

32.0126 

47.1634 

0.6277 

4.7148 

27.0366 

42.0854 

0.6006 

5.2319 

0.9532 

1.8408 

1.2752 

5.2336 

0.9613 

1.6649 

1.3010 

5.2324 

0.9742 

1.6543 

1.3037 

5.7912 

0.7939 

0.7414 

0.8783 

5.8145 

0.7438 

0.7644 

0.7030  . 

5.7949 

0.7872 

0.7281 

0.8540 

6.0145 

0.7866 

0.6897 

0.7882 

5.9967 

0.7946 

0.6878 

0.7870 

6.0252 

0.7931 

0.6867 

0.7860 

5.7059 

0.6251 

0.7353 

0.7195 

5.7159 

0.6240 

0.7326 

0.7174 

5.6926 

0.6109 

0.7270 

0.6707 

5.2319 

0.9532 

1.8408 

1.2752 

5.6862 

0.6148 

0.7374 

0.6710 

5.6926 

0.6109 

0.7270 

0.6707  . 

6.0252 

0.7931 

0.6867 

0.7860 

6.0276 

0.7961 

0.6863 

0.7857 

6.0270 

0.7961 

0.6862 

0.7858 

6.0142 

0.8078 

0.6856 

0.7924 

5.9267 

0.9214 

0.6857 

0.7564 

5.9228 

0.9220 

0.6845 

0.7560  . 

5.9197 

0.9218 

0.6853 

0.7565 

5.8233 

0.9331 

0.6939 

0.7206 

5.7243 

0.8906 

0.7027 

0.7403 

..(Case  2) 


(controller  **) 


(Case  2) 


(Case  1) 


Sensitivity,  and  Weighted  Output  C 
Constraints 


•jin' 


lementarv  Sensitivity 


Case  1  and  Case  2  for  the  4th  and  6th  order  controllers  are  almost  identical,  as  there  is  not 
a  substantial  reduction  of  the  infinity  norms  below  one,  as  will  be  shown  in  Figure  A-l. 


4th  order  controller 
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1.0431 
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6th  order  controller 
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1.0345 

A-5 
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0.9541 
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The  following  figures  map  the  boundaries  for  the  fourth  and  eight  order  controllers  for 
two  and  three  H..  constraints.  As  was  shown  in  Chapter  6,  the  third  constraint  (a 

weighted  output  complementary  sensitivity)  was  not  dominant,  and  therefore  we  can 
combine  the  results  of  the  mixed  controller  with  two  H_  constraints  and  the  controller 
with  three  H_  constraints  of  the  same  order.  This  is  shown  in  Figure  A-l  for  the  4th  order 
controller  and  Figure  A-2  for  the  8th  order  controller. 
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Figure  A-2  Boundary  for  the  eight  order  controller  with  two  and  three  H„ 

constraints 
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